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ABSTRACT 


This study deals with the evaluation of the consequences of the 
release of fission products from a nuclear reactor sited in the ocean. 
The ocean is divided into four distinct zones for which current 
velocity profiles and characteristic diffusion parameters are 
established. These zones are: Zone I (approximate depth of 50 meters 
or less) which is characterized by complete vertical mixing and tidal, 
estaurine, and wind induced currents up to one meter per second; 

Zone II (50 to 200 meters deep) which has as a major feature a 
thermocline located at one-half the depth, below which currents up to 
one meter per second may exist; Zone III (200 to 2000 meters) which 
consists of an upper mixed layer and a deep ocean layer separated by 

a thermocline at a maximum depth of 200 meters; Zone IV (greater than 
2000 meters) which consists of four current layers with bottom 
currents as high as 10 centimeters per second. 

A new fission product release inventory for undersea application 
is developed consisting of 100 percent of the soluble (including 
gaseous fission products) and one percent of the insoluble fission 
products or approximately 66 percent of the gross fission product 
inventory. Based on Carter and Okubo's diffusion model and with the 
assumptions of no current variance, zero mean vertical current 
velocity, and depletion of the inventory by radiological decay only, 
equations are presented which describe the physical transport and 


dispersion of the radioisotopes. 





iv 


A computer program, SEADIF, is applicable to a person immersed 
in the water and is used to determine for both contained and un- 
contained systems the exclusion radius (a gamma dose less than 25 
rem in the first four hours), the low population radius (less than 
25 rem for an infinite time), and the gamma dose rate and radio- 
isotope concentrations as a function of time and position. 

Biological uptake by the ecology leads to concentrations of 
certain radioisotopes in man's food chain in excess of maximum 
permissible concentrations. Results of an investigation of an 
instantaneous release from a 10 MW(t) reactor located in Ocean Zone II 
show that radioisotopes may exist in concentrations above permissible 
limits for as long as 2000 hours after release. A similar investi- 
gation of a continuous release from a 10 MW(t) reactor is carried 


out. 
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CHAPTER! 


INTRODUCTION 


A. Purpose 


An eventual step in man's determination to explore the ocean 
depths, whether for scientific, military, or industrial purposes, will 
be the siting of nuclear reactor systems on the ocean floor. This step 
can safely be taken only after a careful study has been made of the 
possible consequences of an accidental release of fission products to 
the ocean environment through a malfunction of the reactor plant. 

The probability of a major accident in a land based reactor 
plant is exceedingly aN oa and it can be assumed that similar design 
and construction practices will lead to equally remote possibilities 
of an accident for a reactor sited in the ocean. As with reactors 
sited on land, the major concern for undersea reactors is the release 
of fission products to the environment. In this case, the environment 
happens to be the sea instead of the atmosphere. 

For reactors sited on land, it is necessary to investigate the 
effects due to the direct radiation from both contained and airborne 
fission products as well as the effect of the inhalation of the air- 
borne radioactive coneuunneses In the seas, the problem is 
considerably different. Because of the shielding effect of the water, 
the direct radiation effects of a passing cloud of fission products 
are more important than the direct radiation effects of the contained 
fission products. The cloud of contamination is diffusing and moving 
in a marine environment rather than an atmospheric environment. 


Inhalation of radioisotopes is less likely. But the extremely complex 





ecological realm that is interposed between the fission product release 
and man is most significant. In an area of fission product release, 
the ecology will, through biological uptake, concentrate radioactive 
nuclides in their systems, and in so doing, introduce the contamination 
into a primary link in man's food chain. The location of a reactor in 
an undersea environment introduces a number of potential hazards that 


must be considered. 


B. Safeguard Considerations 


1. Direct Gamma Radiation. In the event of an accident, 
personnel in the vicinity of the fission product release may be 
exposed to direct gamma radiation. The primary concern is the direct 
gamma dose received from passage of the radioactive cloud which is 
diffusing in the water and being transported by the ocean currents. 
Because of the shielding effect of the water, the dose due to 
contained fission products is negligible. 

2. Ecological Contamination. Introduction of fission products 
into the ocean leads to contamination of the marine biota (flora and 
fauna). Fish may be exposed to the radioactive cloud, absorb and 
concentrate radioisotopes in their bodies, and then move considerable 
distances before being caught for eventual human consumption. Like- 
wise, shellfish, seaweed, and other commercial ocean assets may become 
contaminated beyond limits acceptable for commercial users. 

3. Contamination of Vessels and Fishing Gear. Vessels passing 
through contaminated waters may introduce radioactive contamination 
throughout the ship's cooling and flushing water systems. In addition, 


the hull and other exposed portions of the vessel may become 
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contaminated. Fishing nets and trawling gear may become contaminated 
with subsequent potential exposure of members of the crew. 

4. Beach Contamination. Ocean currents may carry the contami- 
nation many miles from the point of release. Contaminated waters 
could reach beach, island and coastal areas which are used for 
industrial purposes as well as recreational uses such as swimming, 
boating, or sport fishing. 

5. Bottom Deposition. Contamination may settle out or 
precipitate to the bottom. Commercial dredging or mining operations 
could be affected, as well as the contamination of underwater cables 
which may subsequently be handled in the course of maintenance or 
repair. 

The potential hazards of an accidental release are complicated 
by the fact that reactors may be sited in international waters, and 
the political as well as safety implications of any release will have 
to be considered. 

In general, the contamination of vessels and fishing gear may 
be dealt with by post accident control of access to the affected 
waters. For this reason, it is important that current conditions may 
be accurately predicted from previous oceanographical surveys so that 
the extent of contamination may be estimated and all shipping warned 
of areas in which operations may be hazardous. Current monitoring 
will be required after an accident until surveys show that no 
radiological hazards exist. 

Proper site selection would minimize beach contamination. 
Sites near shorelines should be selected only after extensive ocean- 


ographical surveys have indicated proper offshore current conditions 





exist. An effort should be made to select sites so that minimum 
distances to the shoreline would exceed the distance to which 
concentrations in excess of allowable limits would extend. 

The hazards in conjunction with bottom deposition are more 
easily handled. Again the key is accurately charting the course of 
travel of the contamination. Marine products harvested or mined from 
the bottom are not normaily subject to major horizontal movement. 
Through post accident control, access to the marine products in an 
area of the ocean floor suspected of being contaminated could be 
temporarily suspended until samples had conclusively shown that any 
radiological contamination was less than permissible limits. 

The direct gamma radiation hazard to man and contamination of 
the biota remain as the major problems. A reactor system may be 
the power source for life support facilities which implies that 
personnel will be in the vicinity of the reactor. How close may they 
be? How near may other facilities which house personnel be to the 
reactor complex? To what distances would it be necessary for 
personnel to be evacuated? From what area would it be necessary to 
restrict fishing or shipping? These questions have ied to the 
following definitions applicable to ocean sites. 

EXCLUSION RADIUS (SEA) is the minimum distance at which a 
person immersed in the water with a self contained breathing system 
would receive a total gamma radiation dose to the whole body not 
exceeding 25 rem in the first four hours following the accidental 
release of fission products. 

LOW POPULATION RADIUS (SEA) is the minimum distance that an 


individual immersed in the water with a self contained breathing 





apparatus could be located from the reactor such that he would 
receive not more than 25 rem total dose for an infinite period of 
time. 

SAFETY RADIUS (SEA) is the maximum distance from the reactor 
that concentrations of radioisotopes in excess of maximum permissible 
limits in seawater may exist. The safety radius defines an ocean area 
in which fish would be suspect and contamination of vessels, fishing 
gear, and the ocean floor would be possible. 

The choice of the whole body dose limitation of 25 rem referred 
to in the above definitions is based on the same philosophy as used in 
OMCFR 20° and TID yee This corresponds to the National Committee 
on Radiation Protection recommendations which state that a once in a 
lifetime dose of 25 rem need not be included in the determination of 
the radiation exposure status of a person. 

It may be well to consider the differences between the 
preceding definitions and the corresponding definitions for land based 
reactors (see TID 14844). Three such differences are particularly 
apparent. 

(a) The time for the exclusion radius determination is set 
at four hours instead of two hours. This increase is due to the fact 
that it may require a greater amount of time to evacuate personnel 
from the ocean floor than from the vicinity of a land based reactor. 

(b) No mention is made in regards to thyroid dose due to 
iodine. This concern is generated when the radioactivity is inhaled. 
As it has been assumed that all personnel have an air supply system 
independent of the environment and the reactor compartment, no 


possibility exists for inhalation of radioactivity. 
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(c) Only direct gamma radiation from the passing cloud of 
radioactive material in solution is of consequence. The alpha and 
beta radiation is effectively shielded by the water. 

It is the intent of this paper to present, within the frame- 
work of the preceding definitions, guidelines for evaluating the 
consequences of an accidental release of fission products to an 
ocean environment. In the next chapter, the ocean environment is 
examined and four idealized ocean zones are established. Conservative 
models of diffusion and dispersion are developed for these zones. 

In Chapter 3, the methods are presented for determining the maximum 
permissible concentrations of radioisotopes in seawater. In 

Chapter 4, a fission product release inventory is postulated for two 
cases - a contained (continuous fission product release) and an un- 
contained (instantaneous release) reactor system. An appropriate 
diffusion model is selected and together with the assumptions specified 
in Chapter 5, mathematical models are developed to describe the results 
of the transport and diffusion of the radioisotopes. 

Having specified the model and defined the limits, the results 
of calculations of the exclusion radius, low population radius, and 
safety radius for a 10 MW(t) reactor located on the ocean floor are 
presented in Chapter 6. Conclusions based upon these results and 


recommendations for future work are discussed in Chapter /. 





CHAPTER 2 


THE OCEAN ENVIRONMENT 


In order to determine safe distance factors for undersea 
reactors, it is necessary to have an understanding of the nature of the 
environment which will be the vehicle for the diffusion and transpor- 
tation of any radioactive contamination that would be accidentally 


released. 


A. Topography 


The ocean may be divided into several topographical regions. 
The continental shelf accounts for about eight percent of the ocean 
floor and extends from the coast to depths of about 200 meters with 
an average slope of 0° 07' or about 12 meters per ile This 
submarine terrace, which borders the continents, may extend to 
distances greater than seven hundred miles as is the case of the Sahal 
Shelf off Northern Australia or along the coast of Siberia in the 
North Polar Sea, or it may be essentially non-existent as is the case 
south of Miami. It has an average width of 70 nautical miles off the 
Fast Coast of the United States and an average width of 40 nautical 
miles =rayaliat se The continental shelf breaks sharply at about the 
200 meter depth. 

Extending from this break to the ocean depths, at an average 
slope of 4° 17' for the next 2000 meters of descent, is the 
continental Spee Throughout both the continental shelf and the 


continental slope there is a prevalence of submarine canyons, some 





equal to the Grand Canyon in The origin of these broad floored, 
steep sided troughs has not been fully explained, although recent 
evidence indicates that high velocity turbidity currents could cause 
undersea erosion. The Monterrey Canyon off California is 1600 meters 
Eee pee 

The deep ocean floor consists of oceanic rises, the flat abyssal 
plains, and the mid-ocean ridges. The mid-ocean ridges are inter- 
connected to form a long, nearly continuous ocean mountain system that 
winds over the entire globe. It is interesting to note that over 50Z 
of the earth's surface lies between -2240 meters and -5750 meters and 
that the mean depth of the sea is -3800 peieroe’ 

Rising a minimum of 1000 meters above the deep ocean floor are 
numerous steep sloped underwater mountains or ''seamounts.'' Many of 
these at one time rose above the surface of the ocean where they were 
subject to erosion by the waves. Subsequent subsidence beneath the 
ocean surface has resulted in flat-topped seamounts called "guyots." 
Seamounts are most prevalent in the Pacific Ocean area where over 1400 
have thus far been Biggavsced” They occur in all the other oceans as 
well, but are not nearly so numerous. An interesting feature of sea- 
mounts is the surprising absence of sediments from their summits. 
Photographs have shown rock outcroppings with only small amounts of 
sediment. The bare rock surface and ripple marks in the sediment 
indicate strong currents on the order of 16 centimeters per second 
or more. 

Mid-ocean canyons are found in the abyssal plains. These 
canyons are generally shallow (approximately 200 meters deep), steep 


' ; . ; y sche det 
sided, flat floored depressions, one to five miles wide. Distinct 





from the submarine canyons are the deep ocean trenches which occur 
along the base of about one-half of the continental slopes of the 
world. These may be either V-shaped or have relatively narrow, flat 
floors. Table 1 gives the relative magnitude of some of the world's 


major known deep ocean trenches. 


B. Currents 


The oceans cover 70.8% of the earth's surface and contain about 
331 million cubic miles of water, nearly all of which is in constant 
motion. This motion ranges from near zero currents at great ocean 
depths to velocities in excess of 150 centimeters per second (3 knots) 
in the Gulf Stream. A vertical separation of a few meters can find 
water masses moving in opposite directions. The currents which 
dominate the circulation of the water in the oceans may be categorized 
as (1) wind induced surface currents, (2) major ocean currents, and 
(3) deep ocean currents. 


: The wind is the principle 


Wind Induced Surface Currents. 
current inducing force in the surface and near surface waters. Through 
frictional drag at the wind-water interface, the surface water layer 
is put into motion as is each successive layer of water by the layer 
immediately above. Due to the slippage between the layers, the water 
velocity will be greatest at the surface and will decrease exponentially 
with depth. Because of the Coriolis force which results from the spin 
of the earth on its axis, the current direction is deflected to the 
right of the wind direction in the northern hemisphere and to the left 


of the wind direction in the southern hemisphere. 


As a consequence of the earth's rotation, the surface current is 





TABLE 1. Depths of Deep-Sea Trenches (in meters. 


PACIFIC OCEAN 


Marianas Trench (Challenger Deep), western Pacific 
Tonga Trench, southwest Pacific 

Kurile-Kamchatka Trench, northwest Pacific 
Philippines Trench (Mindanao Trench), western Pacific 
Kermadec Trench, southwest Pacific 

Japan Trench (Idzu-Bonin Trench), western Pacific 
North Solomons Trench, western Pacific 

New Hebrides Trench, southwest Pacific 

Yap Trench, western Pacific 

Palau Trench, western Pacific 

Peru-Chile Trench, southeast Pacific 

Aleutian Trench, north Pacific 

Nansei Shoto Trench, western Pacific 


Middle America Trench, eastern Pacific 


ATLANTIC OCEAN 


Puerto Rico Trench, western Atlantic 


South Sandwich Trench, southwest Atlantic 


INDIAN OCEAN 


Sunda Trench (Java Trench) 


10 


11,034 
10,850 
10,543 
10,033 
10,000 
9,803 
9,142 
9,038 
8,602 
Brie 
8,057 
7,672 
7,512 


6,662 


9,202 


8,262 


7,452 
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directed at an angle of 45° with respect to wind directions. With 
increasing depth, this angle increases and at a certain depth the 
current is directed exactly opposite to the surface. This depth is 
known as the "depth of frictional influence."' At this point, the 
current velocity is approximately 1/23 of the value at the surface. 
The depth of frictional influence is the approximate lower boundary of 
the "surface layer" and is located approximately 50 meters below the 
surface in the middle latitudes and at greater depths at locations 
nearer the equator. A schematic vector representation of this current 
distribution is a logarithmic spiral, known as the Ekman spiral, and 
is shown in Figure l. 

In areas where the ocean depth exceeds the depth of frictional 
influence, the effective mass transport of the water in the surface 
layer is approximately at right angles to the wind direction. 

Surface currents due entirely to wind effects rarely exceed 50 
centimeters per second. For an average wind velocity of 10 meters 
per second, the surface current is about 35 centimeters per second. 
Theoretical determinations of the surface current for a known wind 
velocity depend upon an effective eddy viscosity coefficient and 
are therefore very difficult to Reais ve Large scale measurements 
have shown surface currents to be about 1/33 of the wind speed at a 
height of 30 feet above the Rater However, when associated with 
large scale ocean currents, surface currents may reach 100 centimeters 
per second. An example of this is to be found in the equatorial 
regions where strong steady trade winds exist. 

Major Ocean Currents. The large permanent ocean currents are 


all related to the density distribution of the oceans but are 
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Figure 1. The Ekman Spiral. A schematic =p acne of 
wind induced currents. (Sverdrup ) 


maintained by the prevailing afleiee One dominant feature is the 
circular current motion in each of the major oceans - clockwise in the 
northern hemisphere and counterclockwise in the southern hemisphere. 
The center of this circular motion is offset due to the rotation of 


the earth. 


These current systems can at best be described on an average 





ake 


basis only. Seasonal, as well as day-to-day variations in extent, 
velocity, and direction can be expected. Major current systems such 
as the Gulf Stream or Equatorial Currents are actually narrow, swift 
ribbons separated by counter currents and are subject to considerable 
meanderings. Currents in excess of 100 centimeters per second may 
exist in association with these systems. Table 2 gives representative 
current velocities for a few of the major ocean currents. 

In addition to the wind and density differences, these current 
systems are affected by the topography. The Gulf Stream, which is 
perhaps the most intensively studied current system, extends to depths 
of 1600 to 2000 meters along its eastern edge and is guided by the 
continental slope. 

Deep Ocean Currents. These currents are set in motion largely 
through the movements of cold, dense water which originates in the 
higher latitudes. This water sinks beneath the warmer surface waters 
and spreads out over the vast depths of the oceans. The Antarctic 
Intermediate Water and the Antarctic Bottom Water move northward at 
depths of approximately 1000 meters and 4000 meters, respectively. 

The Intermediate Water moves with a velocity of from six to nine 
centimeters per second whereas the average velocity of the Bottom Water 
is on the order of 0.5 to 1.0 centimeters per second. However, in 

some areas, particularly along the Western Atlantic Trough, velocities 
up to 25 centimeters per second for the Bottom Water have been observed. 
Similar sources of cold water exist in the Arctic, particularly in the 
Greenland Sea between Spitzbergen and Greenland. 


An important layer of deep water movement in the Atlantic has 
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Table 2. Relative Magnitude of Major Ocean Currents 


Magnitude Depth 
(cm/sec) | (meters) 





Atlantic Ocean 


Florida Current 30 — 180 600 
Gulf Stream 10) 500 - 2000 
East Greenland Current 25) Mao) 
Norwegian Current 30 
Straights of Gibraltar 100 | 
Equatorial Countercurrent 30 | 
| 
Pacific Ocean 
Kuroshio JD 0—.180 700 - 1400 
California Current 15 - 30 SO0O7— 1000 
Cromwell Current POO sl 0) | 200 - 400 
South Equatorial Current 50 - 65 
North Equatorial Current | ie > : 
Pacific Equatorial Countercurrent | 30 - 60 | 
Tidal Currents 
San Francisco Bay, maximum ebb 230 | 120 
Bungo Straights 600 


Akutan Pass, Aleutians, maximum ebb | 70 - 450 | 


Note: Data for this table obtained from various sources (refs. 6, 7, 
sea © uel il a lays 
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its origins in the Mediterranean Sea. Dense, saltier water moves out 
of the Mediterranean Sea, over the Gibraltar sill, and spreads out at 
an ever increasing depth. Similar intrusions of high salinity water 
probably occur from the Red Sea and Persian Gulf to provide similar 
sources of deep water movement in the Indian Ocean. 

Actual current measurements in the deep ocean are sparse but 
generally range from 0.05 centimeters per second to 18 centimeters 
per second. Deep water currents at the equator range from 3 to 12 
centimeters per seeeucl. Near Bermuda, measurements from 3 to 18 
centimeters per second and higher have been eran iad Sloe 

Bottom currents were until very recently believed to be near 
zero. However, such evidence as ripple marks and scouring effects 
have been noted in photographs taken of the ocean floor at great 
depths. These indicate significant bottom currents. Theoretical 
calculations for the spreading of Antarctic Bottom Water give 
results as high as 12 centimeters per second along the Western 
Aelantic eerieda A rule of thumb applied in the design of deep 
sea moorings assigns a current of 10 centimeters per second at the 


Mectom. > 


C. Vertical Mixing 


In general, vertical mixing is limited to the surface layer. 
Only in rare instances under the most severe storm conditions would 
vertical mixing effects extend below 200 meters. 

Immediately beneath the mixed surface layer is a zone called 
the thermocline in which the sea water temperature decreases 


rapidly resulting in a rapid density increase. The thermocline, or 
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pycnocline as it is sometimes called, separates the mixed surface 
layer from the layers of intermediate and deep waters. Because of the 
strong density gradient, the thermocline acts as an effective barrier 
to vertical exchange. The thermocline is more pronounced and nearer 
to the surface during summer months. With the cooling of the surface 
waters in winter, it is lower and less distinct. As shown in Figure 
2, the thermocline is less distinct at higher latitudes due to a more 


even temperature and density distribution. 
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Figure 3 shows the clearly developed thermocline that exists 


in the tropical and subtropical ocean areas. Figure 3 also shows the 
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characteristic intermediate salinity maximum which occurs in the trop- 
ical regions of converging surface currents. In these regions, sub- 
surface currents carry high salinity subtropical water under the low 
salinity tropical waters. Surface salinity is a minimum near the 
equator. 

Changes in surface temperature and wind conditions can lead 
to an additional density discontinuity in the surface layer. 
Stratification of the surface layer due to this additional thermocline 
varies on a seasonal as well as a daily basis. 

Since the density of sea water increases with depth, vertical 
random motion is impeded. A mass that is brought to a higher level 
will be surrounded by less dense water and will tend to sink back to 
the level from which it came. Conversely, a water mass moving down- 
ward will be surrounded by denser water and will tend to rise. This 
stable stratification greatly reduces vertical turbulence although 
its effect on horizontal turbulence is negligible. Artifical radio- 
activity introduced at the surface moves vertically downward at a 
velocity of about 0.1 centimeter per second. This motion ceases 
abruptly at the chefnocling = 

Vertical mixing may be increased by an ocean phenomenon known 
as “upwelling" which occurs most commonly in areas near the coast 
where the coastline lies to the left of the wind direction, or in 
areas where major currents diverge. In these areas, there may be a 
movement of the deeper waters towards the surface to replace the 
surface waters which move away from the coastline or diverge because 
of the current patterns. However, it is unlikely that these waters 


would come from depths in excess of 200 to 300 meters. This phenomenon 
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is most conspicuous along the coasts of northwest and southwest Africa, 
California, and Peru where the prevailing winds blow parallel to the 
coasts and because of the Coriolis effect, the lighter surface waters 
are blown aCe Upwelling due to major current divergence is 


most common in equatorial regions. 


D. Site Categorization of the Oceans 


In order to quantitatively study the dispersion of fission 
products in this complex environment, the ocean has been divided into 
four, distinct, depth dependent zones for which current velocity 
profiles have been established and characteristic diffusion parameters 
may be determined. These idealized zones have particular value in 
early site selection surveys when detailed oceanographical data for 
a specific site is not available. They are generalized and must be 
treated as such with every effort made to gain further information 
from site investigation to adapt these zonal descriptions to more 
nearly describe actual conditions. This is particularly true in 
Zone I which includes the shallower ocean areas where local conditions 
may vary a great deal. The velocity profiles associated with these 
zones represent realistic maximums. Average velocity conditions are 
as much as an order of magnitude lower. 

Ocean Zone I (Depths on the order of 50 meters or less) 

Although this zone is characterized by extremely complex 
current patterns which may be influenced by tidal movements and 
estuarine circulation, a uniform current velocity distribution on 
the order of 100 centimeters per second has been established as shown 


in Figure 4. Because the currents are dominately wind induced and 
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the zone is subject to coastal upwelling, it is assumed that complete 
vertical mixing exists throughout the entire depth. 

In these shallow depths, the temperature distribution throughout 
is more variable than is normally expected in the ocean. A thermocline 
may exist at some depth but it would be subject to seasonal, if not 
daily fluctuations. In the winter, the thermocline would be less 
distinct. In any event, it is assumed that it is not an effective 
barrier to vertical mixing. 

At these shallow depths, light penetrates to the bottom to a 
sufficient degree required for rooted plant growth. Because of the 
abundant plant life throughout this depth, marine animals are 
plentiful. This zone has the largest benthic population (flora and 
fauna of the ocean floor) of the entire ne: And together with 
the fact that these areas are generally the easiest ocean areas for 
man to reach, both sport and commercial fishing are to be expected. 
This zone covers approximately four percent of the ocean floor. 

Certain areas fall into Zone I which because of local 
conditions warrant special consideration. This includes consider- 
ation of tidal fluctuations and estuarine conditions. Submarine 
and coastal topography may define inland seas, sounds, gulfs, basins, 
estuaries, harbors, fiords, straights, and other features for which 
particular current patterns and descriptive oceanographical para- 
meters may apply. 

In fact, because of the complex nature of the current patterns 
and the boundary effect of the submarine topography, application of 
Zone I for siting purposes for nuclear reactors must be used 


extremely judiciously. In particular, because of the vertical mixing 
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effect and the shallow depths, the effect of the surface boundary may 
have to be taken into account when applying the equations for 
determination of exclusion radius, low population radius, and the 
safety radius. 

Ocean Zone II (50 to 200 meters deep). Zone II (Figure 5) 
corresponds generally to the continental shelf and although it consists 
of only four percent of the ocean floor, it is a region of high 
interest at this time. Because of the relative shallow depth of 
this zone, enough sunlight penetrates for floating plant life to 
exist at all levels so that marine life is abundant. However, it is 
unusual for rooted plant growth to exist even in the shallower 
regions of this zone. As is the case in Zone I, this area is easily 
accessible and is one of the most widely fished ocean areas. 

The prominent feature of Zone II is the thermocline located at 
one half the depth which imposes an effective barrier to vertical 
mixing. Above the thermocline there is complete vertical mixing 
and high surface currents. These currents are wind induced. Beneath 
the thermocline currents are generally much less, but data is indeed 
meager. Due to the possibility of the influence of major current 
systems, a uniform current velocity distribution of 100 centimeters 
per second has been established. Currents on the order of 15 centi- 
meters per second are more probable. An important point to remember, 
however, is that the direction of the current beneath the thermocline 
may be at any orientation to the surface currents (including opposite). 

Vertical movement beneath the thermocline is assumed to be 
limited to diffusion and transfer by the ecology only since the mean 


vertical current 1S assumed EO be Zero. 
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Ocean Zone III (200 to 2000 meters deep). Zone III corresponds 
to the upper portion of the continental slope and comprises approxi- 
mately 10 percent of the ocean area. [It consists of an upper mixed 
layer and a deep ocean layer separated by a distinct thermocline 
at a maximum depth of 200 meters. The upper mixed layer is similar 
to the corresponding layer in Zone II. It is characterized by 
complete vertical mixing and wind induced currents. There is 
essentially no penetration of the lower boundary (thermocline) except 
by the marine life. 

Beneath the mixed layer, the maximum currents may be as high as 
150 centimeters per second in the area of major ocean currents. This 
maximum velocity will occur very near the lower boundary of the thermo- 
cline. The current will diminish with increasing depth until a zero 
velocity point is reached a short distance from the bottom. Beneath 
this, bottom currents on the order of 15 centimeters per second or 
less flow in the opposite direction. In this region, vertical 
movement is limited to diffusion. Zone III is shown in Figure 6. 

Ocean Zone IV (Depths greater than 2000 meters). Zone IV 
includes about 82 percent of the ocean floor and may be divided into 
four layers as shown in Figure /7. 

1. The mixed surface layer which extends to 200 meters. 

2. The major current strata extending from 200 meters to 
1600 meters deep. 

3. The deep current strata extending from 1600 meters to 
2/3 the total depth. 

4. The bottom current layer which extends to the ocean floor. 


Layers (1) and (2) correspond closely to the layers of Zone III. 








- = = —= aa 
ep» a > 
— = = — i — i —_ 
po ~@ . > >< 
— - ——— so - = 
Gee _ “a a6 
> «¢ =. —- —- 
———/> —- ©£ —» ——— — 


7 = 6 Oe a 





25 


(sz19708WU Q00Z 92 00Z) III su0Z uFZD0 °Q 9ANBTY 


oes/Wwo CT = °A Das/WO OCT = “A “xeW 


VAL MD YL ll A. 






_-— uot ngtz3sT¢ 
querzaing 


PezTTeopy 


nl 


uoTanqTazystqd —-—~ 
quazaing TenjoVy 


a 





srdiew oooz = a Suadeq °xeW 


AZoTooq Aq ueyL 79Y10 
uoTtIeAID.US_g ONES 


eu Ww OOZ 19 Z 


d3s /Wd OOT 





26 


| 
| 
be 
| 
| 
on 
T 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 


(S1910U QOOZ UeYyI 12978918) AT BUOZ uFZDO "yf vin3Ty 





99Ss/WD OT = °A “xeW 


SINATHNO WOLLOd 


deSs/WD CT 


~~ 
_t bee f.———— 


AB OTOOY 
pue uorsnjg std 
Aq uotierjseueg 


SLNHYHNO dead 


0007 49 g 


a) 
e 
z 


°“XebuU 





des/Wo OOT | 98S/wo OCT = °A °xeW 


q ‘yjdeq 








27 


The idealized current distribution in layer (2) runs uniformly from 
150 centimeters per second at the thermocline to 15 centimeters per 
second at the 1600 meter mark. 

The deep currents are represented by a uniform velocity 
distribution of 15 centimeters per second maximum, whereas the bottom 
currents are somewhat lower, with a maximum of 10 centimeters per 
second. Average deep ocean currents applicable to these layers are 


probably as much as an order of magnitude lower. 
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CHAPTER 3 


DETERMINATION OF MAXIMUM PERMISSIBLE CONCENTRATION IN SEAWATER 


One of the potential hazards most difficult to analyze is the 
contamination of the marine life. At present, the upper layers of 
the seas are the source of most fish taken by commercial fisheries 
and eventually consumed by humans. However, some species of food 
fish, although caught in the upper layers, are also found at great 
ocean depths. The commercially valuable sablefish has been found 
in waters 2000 meters deep. Other commercial organisms such as 
deep sea sole, cod, hake, squid, and shrimp are known to live or 
migrate into very deep waters at certain times in their life pelea. 
Centers of fish populations undergo daily as well as seasonal vertical 
migrations. Many species live in water at depths of up to 2000 meters 
during the daytime but at night rise to within 200 meters of the 
surface. These migrants, perhaps only of limited commercial value in 
themselves, constitute a large part of the diet of many food fishes. 
At the present time, fishing techniques and commercial markets 
make fishing at great ocean depths uneconomical. However, one species 
of red shrimp was commercially fished in 400 meters of water off 
Florida in the late tosote Edible varieties of fish exist in 
abundance throughout the ocean depths and await only improved fishing 
techniques to be exploited for commercial consumption. This point is 
especially important as the world's population and food needs continue 


; : 8 
to increase. Acre for acre, the oceans are as fertile as the land, 
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and the potential for development of huge food resources from the sea 
is vast. 

Although human beings do not drink seawater directly in any 
great amount, they do consume many billions of pounds of seafood 
annually. For this reason, if no other, it is necessary to establish 
maximum permissible limits for radioisotopes in seawater. This 
determination of maximum permissible concentration of radioisotopes 
in seawater (MPCC) is complicated by the fact that the biota can 
concentrate certain isotopes within themselves to many times the 
concentration in the surrounding seawater. MPCC's have been 
determined by either the specific activity method or the concentration 
factor method depending upon the degree to which the particular iso- 
tope is absorbed by the body. Before proceeding several terms will be 
defined. 

SPECIFIC ACTIVITY is the ratio of the amount of the radio- 
isotope present to the total amount of the element present. 

CARRIER is a stable element, naturally occuring in seawater, 
which is chemically and biochemically similar to the radioisotope in 
question. 

CONCENTRATION FACTOR is the concentration of the radioisotope 
in the edible portions of seafood divided by the concentration of the 


radioisotope in the seawater. 


meeopecific Activity Method 


For those radioisotopes that are absorbed by the body and for 
which the critical organ is listed by ICRP Publication oe as other 


than the gastrointestinal tract, the specific activity method has 
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been used. This method is based upon the following assumptions: 

1. The radioisotope is introduced into the sea in a form 
such that it acts the same as the stable element or carrier. 

2. The body (both human and marine) is not an isotope 
separator, that is, uptake, distribution, and biological elimination 
is the same for the radioisotope as it is for the stable form of the 
same element. 

3. The chemical makeup of a person who exists on a diet 
consisting substantially of seafood is no different from that of a 
person on a normal diet. 

The MPCC may then be defined as the concentration of the 
radioisotope in seawater such that the specific activity of that 
radioisotope in the ocean is equal to the allowable specific activity 
of that radioisotope in the critical organ. 

The allowable specific activity in the critical organ is 
equal to the maximum permissible burden, q, divided by the amount of 


stable element or carrier in the organ, I The allowable specific 


nb * 
activity in the seawater is equal to the MPCC divided by the 


concentration of the naturally occuring element or carrier in sea- 


icc = SS SS (SL) 


where S = factor to relate the civilian population to the industrial 
workers = 0.1 
X = physical decay constant of the radioisotope 
B = biological decay constant of the carrier 
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Equation (3.1) includes a decay factor which takes into account 
the biological half life of the isotope in man. Because of the 
uncertainties involved in a decay factor for marine species in the 
food chain, such factors were conservatively disregarded. Derivation 
of Equation (3.1) and a more complete discussion will be found in NRC 
Publication 385 and Arcuni's "Survey of Underwater Power Reactor 
Memibicationa."~° Results of the calculations of MPCC's by the specific 


activity method are given in Table C-l of Appendix C. 


B. Concentration Factor Method 

This method is used for determining the MPCC for all isotopes 
for which intake is limited by irradiation of the gastrointestinal 
tract. It is based upon the following assumptions: 

1. An individual derives his entire protein diet from marine 
products. This implies that an individual would eat about 1.5 
kilograms of seafood per week. 

2. The species of seafood consumed is that species having the 
highest concentration factor for a specific radioisotope. 

3. An individual consumes 15 liters of water per week. 

4. The species of seafood remain in an area of contamination 
a sufficient amount of time to come into biological equilibrium with 
its surroundings. 

5. Radioactive decay of the isotope from the time that the 
seafood is caught to the time of consumption is ignored. 

The gastrointestinal tract as the critical organ means that 
the intake is limited by the rate at which radioactive material is 


consumed. The permissible weekly intake is equal to the MPC for 
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water times the amount of water consumed in one week. The maximum 
permissible concentration in seafood is then taken to be this amount 
divided by 1.5 kilograms. The MPCC will be the maximum permissible 


concentration in seafood divided by the concentration factor, f, or 


S x MPC(water) x 15 liters 


ia 1.5 liters’ £ 
MPCC = - (i502) 


where S is equal to the population factor, taken as 1/10 and one 
kilogram of seafood has been assumed equal to one liter of seafood. 

Results of calculations for radioisotopes in which the uptake 
is limited by irradiation of the gastrointestinal tract are given 
in Table C-2 of Appendix C. 

In both methods, it has been assumed that the only source of 
radioactive intake is through marine food products. Any additional 
sources would have to be considered by making appropriate allowances 
for partial permissible concentrations. 

Since, in reality, only a very small segment of the population 
will exist solely on a protein source of seafood, a population factor 
of 1/10 was selected. This is in keeping with the recommendations of 
NBS Handbook 69°. If it were felt that a larger portion of the 


general population could be affected, a factor of 1/100 might be used. 
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CHAPTER 4 


DETERMINATION OF THE FISSION PRODUCT RELEASE INVENTORY 


A maximum hypothetical reactor accident (MHA) is postulated by 
TID 14844 and this accident would be equally applicable to undersea 
locations with only minor - but significant - modifications. Briefly, 
the MHA described in TID 14844 involves Satie complete loss 
of coolant upon complete rupture of a major pipe, with consequent 
expansion of the coolant as flashing steam, meltdown of the fuel 
and partial release of the fission product inventory to the atmos- 
phere of the Mervelscuivanu, 

A fixed reactor system located on or slightly above the ocean 
floor on a supporting structure could be designed as either a 
contained or an uncontained system. For a contained system, the 
envelope surrounding the reactor compartment is designed so as to 
maintain its integrity in the event of a MHA. The envelope of an 
uncontained system is not so designed. 

For an uncontained system, it is assumed that occurrence of 
a MHA results in the inundation of the system by sea water with a 
consequent direct release of fission products to the ocean environ- 
ment. The immersion of the reactor components in seawater is 
important as certain fission products are soluble in or react with 
water. Thus, these soluble fission products are introduced into 
the ocean environment. No reduction in the release inventory is 
assumed to take place due to plating out of fission products on 
reactor component surfaces. Rather, it is assumed that moving masses 


of seawater sweep over the reactor system carrying away the entire 
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fission product release inventory in a single instantaneous "burst" 


type release. These are regarded as conservative model postulates. 

For the contained system, it is assumed that the hydrostatic 
pressure is greater than the containment pressure. During normal 
operation, this means a certain amount of leakage of sea water into 
the reactor structure - a perfectly normal situation which would be 
routinely handled by pumps. After the MHA, even though the con- 
tainment vessel remained intact, it is assumed that power to these 
pumps is interrupted and eventual flooding of the containment 
structure results. The reactor components are therefore immersed in 
seawater and the soluble fission products are released to the flooded 
containment structure. Because the containment structure cannot be 
absolutely leaktight, it is assumed that the fission products 
released to the containment leak out to the environment at a con- 
stant rate of 0.1 percent per day. Barring a major collapse of the 
containment structure, this figure would be extremely conservative. 

The fundamental assumptions which form the basis for 
determining the fission product release inventory are given below. 

CASE I. UNCONTAINED SYSTEM 

1. It is assumed that complete inundation of the reactor 

system by seawater occurs at the time of the MHA. 

2. As a result of this inundation, 100% of the soluble 

fission products and one percent of the insoluble fission 

products are released to the ocean environment. 

3. No plating out or adhesion to the reactor component 


surfaces by the soluble fission products occurs. 
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4. Because of their high solubility in cold water, all 

fission product gases are included in the release inventory 

and it is assumed that none escape to the ocean surface. 

5. The fission product release inventory mixture is assumed to 
decay at the same rate as the gross fission product inventory 


(see Appendix B, Part 5). 


CASE II. CONTAINED SYSTEM 

The assumptions of Case I apply with the addition of the two 
assumptions given below. 

6. The leakage from the containment structure occurs at the 

constant rate of 0.1% per day. 

7. Physical integrity of the containment structure is main- 

tained throughout the time of interest. The time of interest 

extends from the time of the MHA until such time as the 

radioactive decay has reduced the fission product release 


inventory to negligible levels of activity. 


The differences in resultant fission product release inventories 
for land based and undersea reactors are highlighted by Table 3. The 
assumption that 100% of the solubles are released to the environment 
is extremely conservative and accounts for the significant increase 
in the undersea reactor release inventory. However, this conservative 
assumption is felt to be the only valid approach until future work in 
this area provides a sound basis for modification. 

The assumption that 100% of the gases are dissolved in the 


seawater is conservative in so far as the effects on a person immersed 
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in the water is concerned. This assumption is not conservative if the 
investigation is directed to the effects of airborne radioactivity on 
a person on or above the surface of the ocean. Airborne radioactivity 
may result from a gaseous fission product bubble rising to the surface. 
This is more of a possibility in shallow ocean areas and under in- 
stantaneous release conditions. A detailed analysis of the bubble 


effect or rising plume effect was not performed. 


Table 3. Comparison of Fission Product Release Inventories 
for Contained Reactor Systems? 






PERCENT RELEASED PERCENT RELEASED 
TO CONTAINMENT ae. TO ENVIRONMENT 


TID 14844| UNDERSEA j MILD, “een UNDERSEA 


Noble Gases 100 
Halogens 100 
Other Soluble 100 
Solids 
Insoluble Solids 1 af if 
Gross Fission 
Product Invento 66 1B 66 
a Percent figures based on total curies (beta and gamma emitters 
included) (see Appendix D). 
b If the elemental state is soluble in or decomposes in water, then 


the total activity due to that fission product is included in 
the category of soluble. 
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CHAPTER 5 


TRANSPORT AND DISPERSION OF FISSION PRODUCTS 
IN AN OCEAN ENVIRONMENT 


A. Physical Processes of Transport and Dispersion 


After the fission product release inventory has been intro- 
duced into the ocean environment, it will be subject to several 
physical processes. These are (1) molecular diffusion; (2) turbulent 
diffusion; (3) physical displacement by large scale currents; (4) 
shear effect; (5) depletion through precipitation, adhesion to the 
ocean bottom, or uptake by the ecology, and (6) radiological decay. 

Molecular Diffusion. Molecular diffusion is due to the random 
motion of the molecules of liquid and dissolved substances. It is a 
function of the temperature, pressure, and concentration gradients. 
However, in comparison with the effects of the turbulent diffusion, 
molecular diffusion appears negligible. Although it is generally 
neglected in all but extraordinarily calm eairiene.’, orube 
points out that the interaction between molecular and turbulent 
diffusion may be important even in large scale phenomenon. 

Turbulent Diffusion. Turbulent diffusion arises from the 
large and small scale eddies which may be considered as randomly 
distributed in both time and space. The eddies span many scales 
in size and their effect is perhaps best described by Neumann. ''For 
a given volume filled with material to be diffused, certain scales of 
motion serve only to transport the whole volume, others move the 


material around within the volume without spreading it further, and 
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only a certain intermediate range is effective in spreading the 
elements of the volume further peregno 

One of the complicating factors when considering turbulent 
diffusion in an ocean environment arises from the fact that the 
horizontal coefficient of eddy diffusion is proportional to the scale 
of the disturbance. Weigel illustrates this by pointing out that 
the rate at which two particles tend to separate increases with 
increasing distance. At close distances, these particles are moved 
apart by small scale turbulence only, whereas at greater distances 
their separation is effected by the larger ies.” 

For continuous release of fission products from a point source, 
the distribution at any particular instant could be represented by a 
thin ribbon of concentrated fission products meandering through a 
segment of a circle as shown in Figure 8. At a point downstream from 
the point of release, the average concentration measured over a 
period of time would be lower than the concentration in the ribbon 
because of this meandering. This meandering occurs on all scales of 
oceanic eddies, and it is important to remember that it is the 
average distribution rather than the instantaneous conditions which 
should be observed when concerned with the total radiation dose 
received. The gradual widening of the narrow ribbons is due to the 
smaller isotropic eddies, whereas the meandering of the ribbon is due 
to the intermediate range eddies. The general orientation of the 
ribbon is, of course, due to the large eddies (mean current). 

Physical Transport. Eddies of a scale considerably larger 
than the area of interest may be considered as anisotropic entities 


and treated individually. These large scale eddies which are found in 
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of Dye from a Continuous Source. (Okubol®) 


the region of interest will tend to move the escaping fission products 
in accordance with their mean velocity components. The ocean zones 
established in Chapter 2 allow quantitative treatment of these 
characteristic currents. 

Shear Effect. Because the large scale eddies are anisotropic, 
that is, the characteristic length scale is longer in the direction of 
mean flow, a nonuniform velocity profile may appear. This velocity 
gradient combined with vertical diffusion leads to an effective 


SEDO This, in turn, 


horizontal dispersion known as shear effect. 
leads to elongation of the diffusing cloud of fission products along 


the axis parallel to the direction of mean transport. 
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One of the results of poedeaiee investigation of shear effect 
was that the effective coefficient of horizontal diffusion is inversely 
proportional to the coefficient of vertical eddy diffusion. That is 
to say, a stable density gradient increases the horizontal mixing 
considerably. 

Depletion. In addition to radiological decay, reduction of 
the overall fission product release inventory may occur through any 
of the following processes: 

1. chemical precipitation, 

2. absorption to the bottom sediments, and 

3. biological uptake by the biota. 

Chemical precipitation may occur as a result of the inter- 
action of certain fission products with the naturally occuring 
substances in seawater. The resultant insoluble may then settle to 
the ocean floor. The process of co-precipitation, whereby similar 
elements are simultaneously removed from solution, may also occur. 
For instance, strontium is usually co-precipitated and carried along 
with calcium carbonate. 

An even greater portion of the release inventory may be 
removed from the release cloud through sorptive processes by the 
bottom sediments. Finely divided organic and inorganic solids may be 
held in suspension in a layer a few centimeters thick at the ocean 
floor. These solids have a very definite tendency to absorb (remove 
from solution) and to hold substances dissolved in seawater. 

Uptake of fission products by the biota may take place through 
adsorption to the surface areas, engulfment, or through metabolic 


1 e @ e e J 
processes. The ability of certain marine organisms to concentrate 
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radionuclides was discussed in Chapter 3. Subsequent movements by 

the biota can also contribute to horizontal and vertical dispersion 

of the inventory. This is particularly true in the case of the 

extreme vertical migrations of zooplankton which occurs daily. 
Radiological Decay. Radiological decay is the only process for 

which a quantitative estimate of the depletion of the release inven- 


tory may be made with any degree of reliability. 


B. Fundamental Assumptions 


The fundamental assumptions which form the basis for deter- 
mining the extent of the transportation and diffusion of the fission 
product release inventory are as follows: 

Diffusion Model. Diffusion of the fission product release 
inventory is assumed to occur in accordance with the three di- 
mensional diffusion model developed by Carter and Oleupewier 
Recommended by eum, this model includes shear effect and most 
accurately predicts the diffusion process over the entire time 
scale. This model includes the effect of turbulent diffusion due to 
the large scale eddies in both the horizontal and vertical planes. 
For completeness, the assumptions used to develop this model include 
the following: 

a. The horizontal and vertical shear is assumed constant. 

b. The eddy diffusion coefficients in the x-, y-, and z- 

directions are assumed constant. 

c. The initial concentration is assumed to be a point source. 

The first two assumptions do not describe the actual conditions 


in a turbulent ocean. The actual velocity profile is not linear as 
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would be required for constant shear terms. Likewise, in turbulent 
transport, the unsteady motions and eddies of varying scale that are 
doing the transporting are by no means constant. However, at best, 
only an average diffusion coefficient can logically be used to define 
this motion. The basic concept of a diffusion model incorporates 
average parameters to describe the results of random, turbulent 
movement. 

The assumption of a point source does not accurately describe 
conditions of leakage from a containment structure. However, this 
shortcoming may be corrected by using a virtual distance concept. 
This has been done and results are given in Appendix B, Part 10. 

The model, given by Equation (Bl.3) of Appendix B, represents 
the contours of concentration as a set of elongated ellipsoids with 
common principle axes. The major axis is aligned in the direction of 
the mean flow. The lengths of these axes vary with time. The degree 
of elongation depends, apart from time, on the shears, and ae and 
the eddy diffusivities, A. 7 and Ao: The combined effect of the 
shear and diffusion accelerates the rate of disperson of dye around 
the center of mass, i.e., the peak concentration in the mode 

Current Variance. It is assumed that the velocity and direction 
of the mean current does not shift during the time of interest. If 
such shifts occured, a reduction in the concentration could be 
expected. 

Vertical Velocity Component. The mean velocity component in the 
vertical direction is assumed to be zero. Vertical mixing by the 
biota and vertical mixing due to thermal effects is neglected. It is 


assumed that stable stratification exists and vertical mixing is due 
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to the processes of diffusion only. Any vertical mixing that did 
occur would reduce the concentration. 

Inventory Depletion. Depletion of the fission product release 
inventory by precipitation, absorption to the bottom, or uptake by 
the biota is neglected. Any actual depletion that occurs would lead 
to an overall reduction of the concentration. However, particularly 
in the shallow depths of Zone I, concentration of fission products 
on the bottom by these processes could lead to conditions of greater 
importance. In these ocean floor areas which may be inhabited by 
commercial marine products, such as oysters, clams, mussels, and 
other benthic organisms, the equations presented in Section C of this 
chapter may not be applicable. 

Fission Product Decay. A generalized model for the decay of 
fission products is assumed. In this model, no decay is assumed 
during the first 1000 seconds after shutdown. After this time, it is 
assumed that the fission product release inventory decays as k eae 
The assumption of no decay during the first 1000 seconds after shut- 
down may at first appear too conservative. However, it is in part 
offset by the fact that fission products with a half life of less 
than four minutes are not included in the release inventory. The 
overall effect of these assumptions is conservative in that the 
decay rate for the higher energies of the gamma ray spectrum fall off 
faster than do the lower energies with a resultant softening of the 
spectrum with time. 

Reflection from the Ocean Floor. It is assumed that the 
reflective effect from the ocean floor can be treated by adding the 


image with respect to the ocean floor. This is conservative in that 
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adhesion by the bottom sediments would result in considerably less 
than 100% reflection. 

Boundary Effects. The boundary effect of the sea surface and 
the thermocline is neglected. For Ocean Zones II, III, and IV, such 
an assumption is reasonable due to the extremely limited vertical 
diffusion. However, for Zone I, this assumption is not conservative. 
Rather the reflection at the thermocline or the sea surface would lead 
to higher concentrations. 

Steady State Conditions. In the continuous release case, it is 
assumed that steady state conditions are established instantaneously. 
This is extremely conservative in the investigation of locations a 
great distance from the point of release with slow currents. For 
points near the release or for high currents, the assumption is 
reasonable, and it greatly simplifies the mathematics. 

Shielding Effect. In determining direct gamma dose, the effect 
due to the contained fission products is neglected and only the gamma 
dose due to the cloud passage is considered. Because of the shielding 
effect of the water, at distances greater than 10 meters from the 
containment, the contribution of the fission products within the 
containment towards the total dose received is negligible. Also, 
because of the shielding effect of the water and the fact that 
inhalation is not possible, gamma radiation of less than 0.1 Mev 
energy, all alpha, and all beta radiation is neglected in dose 


determinations. 


C. Mathematical Models 


With the assumptions as presented in the preceding section, it 
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is possible to derive equations which describe the physical transport 
and dispersion of the radioisotopes that are assumed to be released 
in the event of a MHA at a reactor sited on the ocean floor. This 
has been done in Appendix B for both the contained and uncontained 
reactor systems. 

For the uncontained reactor system, the total (integrated) dose 
received as a function of distance from the reactor and mean current 


velocity is given by the following equation: 


10° sec + 
' -1.2 
Exo. ve) = DCrexeevsy dt + kt DiGCEs scene) GCE (5.19 
j 
0 Ge 
where D(t,x.,v.) = dose rate due to fission products in the water as 


a function of time after release, distance from 
the reactor, and mean current velocity. 


time after release 


lie 
It 


x 
It 


normalization constant for the decay of mixed fission 
products 


This equation forms the basis for determining the exclusion radius and 
the low population radius - the only difference being the selection of 
the upper limit, T, as four hours or as infinity. The exclusion radius 
or low population radius is the value of x, such that the dose evalu- 
ated by Equation (5.1) is 25 rem. The equation which forms the basis 
for determining the safety radius is given by the following 


expression: 


MPCC = —————_____— (522) 
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where 
1 
C= Saat (3) 
gro! 2 ( Aaa yl/2 
xX 2 
A A 
sa Oa aaa 2 2 
g = 1D cm A + 2 A ) Ce) 
x x 
and M. = the total activity (curies) of radioisotope i, released 
to the environment. 
r. = decay constant for radioisotope i. 


Equation (5.2) may be solved for t, the time after release at which 
the peak concentration of the isotope equals the MPCC. By knowing 
the mean current velocity, a safety radius may then be established. 
For the contained reactor system, the applicable equation for 
exclusion radius or low population radius calculations is given by 


Equation (5.5) where T equals 4 hours, or T equals infinity. 


dt 
3 “1.25, 
ics. V5) seam 10° Sseepaep" (x, Vv.) + D' Cxs,v5) k t (525) 
0 
where D'(x,,v.) = instantaneous steady state dose rate as a function 
of distance from the reactor and mean current 
velocity. 


As in the uncontained case, the exclusion radius or low population 
radius is the value of x, such that the dose, when evaluated by 
Equation (5.5) is equal to 25 rem. 

For safety radius determinations, Equation (5.6) applies. 


co 


MPCC = M, Le W (t,X0,Vo)dt (5.6) 
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where Q = leak rate 


function defined by Equation (B6.1) which 
is the basic mathematical diffusion model 
as determined by Carter and Okubo (Ref. 27). 


WiCEseceniie) 


Because of the instantaneous, steady state assumption, the 
highest concentrations (and therefore the limiting concentrations) 
occur at t = 0, or before radioactive decay begins. The safety radius 
is that value of x, such that the concentration as evaluated by 


Equation (5.6) at t = 0 is equal to the MPCC. 
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CHAPTER 6 


RESULTS 


A computer program, SEADIF, has been written which wiil 
calculate the exclusion radius, low population radius, the direct 
gamma dose rate four hours after release, and the time (or in the 
contained case the distance) required for the peak concentration of 
each fission product radioisotope to fall below MPCC. SEADIF consists 
of two separate and complete programs, the use of which is determined 
by whether it is a contained or uncontained reactor system. The 
following data must be read in by the computer in proper format: 

a. Power level of the reactor 

b. Diffusion and shear coefficients 

c. Current velocity 

d. Distance from the containment 

e. Containment diameter (optional input) 

f. Isotope information for all fission products of 

consequence, including amount per megawatt, MPCC, 
half life, percent decay by gamma energy group, 
and solubility information 

g. Containment leakage rate (contained case only). 


The computer program and input data format is present in Appendix E. 


A. Uncontained System 


The exclusion radius, by definition, depends upon the total 
dose received in the first four hours after initial release. This 


total dose will be a function of the distance away and the mean 
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current velocity. Figure 9 depicts the four hour dose as computed by 
SEADIF for typical current velocities and diffusion coefficients. There 
are several interesting features of this set of curves. 

First, for high velocities and relatively short distances from 
the release point, the total dose received is independent of the mean 
current velocity and inversely proportional to the distance. Although 
higher velocity conditions will mean that a higher concent ration of 
fission products will be carried a certain distance, this concentration 
passes over the point of interest in a shorter time with the net result 
being an identical total dose received. 

Velocity dependence enters the picture when the fission product 
cloud (defined as extending a distance of 30) has not completely passed 
over the point in question within 1000 seconds after release. This is 
because at t equal to 1000 seconds, it has been assumed that radio- 
active decay begins. It now depends upon the time at which the cloud 
passes over as to the total dose received. 

A second interesting feature of Figure 9 is the departure from 
a linear representation on a log-log plot. This occurs when the 
fission product cloud has not completely passed over the point in 
question after four hours. 

If radiological decay is ignored, the results will be as shown 
in Figure 10. Comparison of Figures 9 and 10 illustrate the effect of 
radiological decay. As is to be expected because of the basic decay 
assumptions, if the entire cloud passes over the point in question 
within 1000 seconds after release, there is no difference in the dose 
received. Figures 9 and 10 are based on the same fission product 


release inventory, that is, fission products with a half life of 


















= > o << * 
. J cnt need — —a > =-> EE 
> = [as —> > 
a6 ae ’' 
l= => —_> —_ = Ss | = >= Sa 
_ eo a -m-¢< 
——~ es aa a eo—_—_ z=. —_ o 
a ——e & 
== - - = 
Pe a _ Ss f a - ‘i _, ——_> 2 > 
Soe a= = e“- Se = Ce 41 eS 


Ce 
2 EP kage) OM” cam ——— 
-— 6 eee eee ee 6 Le ee ae 
- a 
re a i a ae 


i. iain —— a 
























en ; 















_ 


Four Hour Dose (rem) 


50 


10 








vo = Current Velocity, meters/sec 


EXCLUSION RB 


10 100 1000 10,000 
Distance (meters) 
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less than four minutes are not included in either case. 

One of the major considerations with regard to removal of 
personnel from the area of an undersea reactor will be the dose rate 
"picture'' at four hours, an example of which is shown in Figure ll. 
For a 10 MW(t) reactor and a mean current velocity of 20 centimeters 
per second the exclusion radius is 5/0 meters. From Figure 1l it may 
be seen, however, that at 5/0 meters the dose rate at four hours is 
essentially zero since the entire cloud has passed. Therefore, after 
four hours no additional dose would be received by personnel within 
that distance from the reactor. 

This points up a very important concept. For an uncontained 
reactor system the low population distance is, except for extremely 
small mean current velocities, the same as the exclusion radius. 
Again, this is due to the fact that once the cloud has passed over the 
point in question the dose rate in that vicinity is negligible. This 
concept is different from the concept presented in TID 14844 because 
of the fact that in an atmospheric environment the dose due to the 
fission products in the containment building is significant. 

Because of the shielding effect of the water this is not the case for 
an undersea location. Figure 12 illustrates the difference between 
the four hour dose and the infinite dose, these doses forming the 
basis for the exclusion radius and low population radius respectively. 
It is seen that differences are apparent only at total doses well 
below the 25 rem limit. 

A great deal of uncertainty lies in the selection of appropriate 
diffusion and shear coefficients. Two representative diffusion 


coefficient data sets were chosen based upon parameters presented by 
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aps, 


Carter and olives. These two data sets, shown in Table 4, define a 
range of realistic diffusion and shear coefficients but should not be 
thought of as totally inclusive. Unless otherwise noted, all figures 
are based upon data set 1. The effect of changing data sets upon the 
exclusion radius is shown in Figure 13. One important difference in 
these two data sets is in the selection of the vertical diffusion 
coefficient, A. Figure 14 shows the effect of varying this coeffi- 
cient alone over the range of values presented by Carter and Okubo. 


| 
Table 4. Diffusion Coefficient Data Sets 
| 


Diffusion Coefficient, x-axis, A, 0.39 faeces 0.39 ic yeee 
Diffusion Coefficient, y-axis, fy OAS diane ele Oe se neyeee 
Diffusion Coefficient, z-axis, A, 1 OeNON = no gee Once lam nae 
oe -6 -1 -4 -1 
Shear Coefficient, ie Iter NG sec 1.0x10 sec 
eae -3 -1 -3 -1 
Shear Coefficient, 2, 1.8x10 ~ sec 6.6x10 ~ sec ~ | 





Results are also sensitive to the horizontal diffusion coeffi- 
cient, AL» as shown in Figure 15. This coefficient is much larger 
than A, because of the much greater diffusion in the horizontal 
directions. 

The effect of the shear terms is not nearly so noticeable. 
Shear coefficients enter into the calculations of 8 only (defined by 
Equation (5.4)) and the exclusion radius results are rather insensitive 
to this term. However, where these terms do effect results is in the 


peak concentration calculations. The larger that 2 is, the greater 
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the elongation of the diffusing cloud. This has the effect of 
reducing the peak concentration. 


ye 


For times much less than 8 it may be seen from Equation 


SY 


(B6.9) that the peak concentration goes as t_ whereas for times 


=1/2 : , : 
much greater than 8 / the peak concentration is proportional to 


a Therefore, by varying 8 the rate of change of the peak con- 
centration can change. Beta is affected by both diffusion and shear 
coefficients. Figure 16, which shows the change in peak concentration 
with time, has been normalized so that the peak concentration is equal 
to l at 10° seconds. Figure 16 is applicable only to long lived 
fission products in that radioactive decay has not been included. 

The reactor power enters into the calculations by way of the 
fission product release inventory. The dose rate and total dose 
received are directly proportional to the power. The results for the 
exclusion radius are graphically represented in Figure l/. 

Up to this point, the effect of the surface boundary and the 
extent of vertical diffusion has largely been ignored. The dose rate 
as a function of distance above the ocean floor was determined at 
times of four hours and eight hours after instantaneous release for 
a 10 MW(t) reactor using diffusion coefficient data set 1. Results 
are shown in Figure 18. It may be seen that as long as the sea 
surface or thermocline is greater than about twenty meters above the 
release point the assumption that the boundary effect may be 
neglected is valid for a time frame of a few hours. 

More important, as long as the reactor location is not in the 
vicinity of upwelling or in shallow waters where vertical mixing 


extends to the bottom, it is not likely that short half life 
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radioisotopes, in harmful amounts, will extend to the mixed layer. 
Therefore, contamination of surface vessels by radioisotopes with a 
short half life is not likely. 

For an uncontained 10 MW(t) system, the critical isotope in 
determining the safety radius is ieee as shown in Figure 19. Using 
diffusion coefficient data set 1 (shaded bars on Figure 19) the time 
required for the concentration of feces to fall below MPCC is 2,130 
hours. For a mean current velocity of 1 meter per second the potential 


safety radius is extremely large. Under these conditions, the 


assumption of no current variance becomes extremely conservative. 


B. Contained Systems 

Despite the conservative assumption that steady state 
conditions are achieved instantaneously it may be seen from Figure 20 
that the four hour dose figures for a contained 10 MW(t) system are 
extremely small. Furthermore, these figures are not as velocity 
dependent as in the case of uncontained systems. 

Briefly examining the steady state assumption, it is noted that 
this is most conservative for extremely small velocities and large 
distances. For instance, at 1000 meters and 1 centimeter per second 
current velocity conditions, the initial activity released in the 
course of a MHA would not begin to arrive until approximately 10° 
seconds after initiation of the MHA. Therefore, the four hour dose 
would be zero. However, Figure 20 reflects a four hour dose for 
these conditions in excess of ton roentgen. 

The main point is that before serious consequences due to the 


direct gamma dose are experienced, an increase by a factor of 100 to 





Time After Release, hours 


Safety Ra 


ae ee —— $s ae oe — aw = 


3000 


Shaded: 


Unshaded: 


dius Determined by te™ 


Data Set l 


Data Set 2 


Safety Radius 


64 


a” a oo aa eo a, ne ee 


2000 


1000 
Np? 


Figure 19. 


3m 


Determined by Te 


iLL il 


RyL?? Rye 0® Tel? /™ Tet?7m Ba 140 Pr 


Critical Isotopes 
Time Required for Concentration of Critical 


Isotopes to Fall Below MPCC for an Uncontained 
10 MW(t) Reactor 


ot44 





Four Hour Dose (rem) 


65 





ror 
—2 
10 ene ORM / SEC 
vo = 0.01 
Analytical Solution ea 
~3 for v, = 0 — 
10 
vo = 0.001 
least 
10 100 1000 10,000 


Distance (including virtual distance), m 


Figure 20. Four Hour Dose as a Function of Distance 
(with decay) for a Contained 10 MW(t) Reactor 
with Diffusion Coefficient Data Set l 





66 


1000 would have to be incurred. An exclusion radius or low popu- 
lation radius has no real significance for a contained reactor system 
until power levels in excess of 1000 MW(t) or leak rates considerably 
in excess of 0.1 percent per day of the fission product release 
inventory are realized. A comparison of the exclusion radius and low 


population radius for a contained and uncontained system is shown in 


iable 5. 


Table 5. Siting Distances (in meters) for 10 MW(t) 
Reactor Located in Ocean Zone II 


Contained Uncontained | TID 14844 





Exclusion Radius 


Low Population Radius 


Low population radius determined on basis of iodine 
dose to thyroid. Limiting distance on basis of 
direct gamma dose is 250 meters. 


Figures 21 and 22 show the effect of the variation of the 
diffusion coefficients on the four hour dose for a contained reactor 
system. 

The principal concern in contained reactor systems is the 
safety radius determination. Results of calculations for a 10 MW(t) 
System are shown in Figure 23. The critical isotope for Ocean Zone II 
with a mean current velocity of 1 meter per second is Nbeaw with a very 
long half life (32,400 hours) and a very low natural abundance of its 
carrier, the element Niobium, in seawater (1 x one grams per kilo- 


gram). These two factors outweigh the fact that the fission product 
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Figure 22, 


Effect of Variation of Horizontal Diffusion 
Coefficient, A,, on 4 Hour Dose for a 
Contained 10 MW(t) Reactor 
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release inventory for eee is only 75 curies per megawatt because of 
its tmselubility. 

An examination of the critical isotope information contained in 
Table 6 points out the difficulty in making any general remarks as to 
what makes an isotope critical. The only common factor among the nine 
most critical isotopes is the relatively long half life of each. 
Other than this, the criticality of an isotope depends upon the MPCC 
and the release inventory. These quantities depend upon one or more 
of a number of factors such as solubility, critical organ, natural 
abundance of the carrier in seawater, concentration factor, fission 
product inventory, and permissible body burden. 

The most critical isotope which determines the safety radius 
varies with release conditions and ocean parameters. For a 10 MW(t) 


reactor, the critical isotope for a contained system with a mean 


pou, 


current velocity of one meter per second is Nb ; for an uncontained 
; ices 129m 
system with the same current profile it is Te ; while for a 


contained system and a current of two meters per second, it is 


Tet /™. 
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CHAPTER 7 


CONCLUSIONS AND RECOMMENDATIONS 


Use of the computer program SEADIF allows a preliminary 
investigation of the suitability of an undersea site for a nuclear 
reactor in so far as radiation hazards are concerned. Results show 
the desirability of containment not so much because of the direct 
gamma radiation hazard but rather because of the possibility that 
exists for contaminating the marine biota in a large ocean area. 
With proper containment, the direct radiation hazard is essentially 
eliminated and the ocean area which could possibly be contaminated 
is minimized. 

The models developed for both the contained and uncontained 
reactor systems are conservative. This conservative approach was 
taken at every step in the development of these models with the end 
result that figures for the exclusion radius, low population radius, 
and safety radius more nearly provide maximum hypothetical limits 
rather than good engineering estimates. A great deal of work remains 
to be done to refine these models. 

Ocean Zone I was defined and then largely avoided because of 
the complex environment and range of conditions that exist. A more 
detailed study of this zone is required with particular emphasis on 
tidal effects, estuarine conditions, shore line contamination, and 
descriptive parameters in order to meaningfully apply the diffusion 
models. This zone will require further subdivision and the effects of 
vertical mixing and surface reflection will have to be considered. 


For all the ocean zones, a need exists for realistic diffusion 
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and shear coefficients. Recent work by the Chesapeake Bay 


VIREDEAS 


Institute is a step in the right direction but much more 
study and experimentation is required. 

Closely allied with diffusion coefficients is the effect of 
current variance. An extremely conservative assumption that was made 
was to assume no current variance. Incorporating the effects of such 
variance, particularly when dealing with the safety radius, would 
result in much more realistic conditions. 

The assumptions involved in the determination of the fission 
product release inventory are based on the best information available 
with the desire for conservatism uppermost in mind. Much needs to be 
done in determining the amount of fission products that would be 
released to the marine environment in the event of a MHA. The 
solubility of fission products in seawater should be investigated, as 
should the chemical and physical processes involved in the removal of 
fission products by precipitation, adhesion, absorption, and uptake 
by the marine life. 

Safety radius determinations have led to the identification of 
a number of critical isotopes. Starting with the isotopes of yttriun, 
niobium, ruthenium, tellurium, cesium, barium, and praseodymium a 
review of the input data used in calculating MPCC's should be made. 
More suitable carriers may be identified and more realistic concen- 
tration factors may be applicable. Marine biological assimilation 
and elimination rates should be investigated. Both the specific 
activity method and the concentration factor method of calculating 


MPCC's have shortcomings. A thorough study in this area could be of 


great benefit. 
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Contamination of the marine environment by activation products 
was not discussed. Corrosion products such as isotopes of iron, zinc, 
cobalt, manganese, chromium, copper, and nickel may be of particular 
interest as these elements have an extremely low natural abundance in 


seawater and high concentration factors in seafood. 
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APPENDIX A 


LIST OF SYMBOLS AND ABBREVIATIONS 


effective cross sectional area of diffusion cloud, 
meters 


eddy diffusivity in the x-direction, ee 
eddy diffusivity in the y-direction, Eo eee 
eddy diffusivity in the z-direction, Pvece 
buildup coefficient 
buildup coefficient 


biological 
man, days 


decay constant of the carrier element for 


dose rate, rem/hour 
ocean depth, meters 
containment diameter, meters 
decomposes in cold water 
energy per gamma photon, Mev 

th 
average gamma energy of the n energy group, Mev 


é © vi 
GrrOr LcuNncEloOn derineda as ——— e 


complimentary error function = 1 - Erf z 


ratio of consumption of water to consumption of 
seafood 


concentration factor of seafood 
total dose received, rem 


total quantity of the carrier in the critical organ, 
grams 





ns 


nw 


insol 


MPCC 


To 


Vo 
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concentration of the carrier in the edible portions 
of seafood, grams/kg 


concentration of carrier in seawater, gram/kg 
insoluble in cold water 


normalization constant for decay of mixed fission 
PLOdUcES «5 70x 102 secl-¢ 


-l 
leak rate, sec 


total activity of mixed fission products released to 
containment, curies 


total activity of isotope i, released to containment, 
curies 


activity of mixed fission products released to contain- 
ment falling into the n-th gamma energy group, curies 


maximum hypothetical accident 


maximum permissible concentration in drinking water for 
168 hour week, uc/kg 


maximum permissible concentration in seawater, uc/kg 


maximum permissible quantity of radioisotope in the 
critical organ, uc 


roentgen 

factor to relate the recommendation for industrial 
workers to the civilian population, equal to = 
concentration, uc/kg 

concentration of a specific isotope i, uc/kg 
peak concentration, uc/kg 

time after release, sec 

time, seconds 

time (dummy variable) 


virtual distance, meters 


mean horizontal current velocity, meters/sec 


2 Ss oe 


we, 





Xo 


ve 


02 
Z, 
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distance measured from peak concentration parallel to 
mean current flow, meters 


distance measured from initial point of release 


distance normal to mean current flow, parallel to 
ocean floor, measured from point of release, meters 


distance above ecean floor (point of release), meters 


3/2 


diffusion term, sec. per meter 
buildup coefficient 
buildup coefficient 


diffusion-shear term, sec 


gamma function defined as p21 = 


gamma energy group 

; : : SE Bll 

incomplete gamma function defined as | t dt 
0 


delta function defined to be zero for all values of x 
except x = 0 


— th 
energy coefficient for n energy group, R-meters 


energy dependent conversion coefficient, R/hr per 
Mev/cm -sec 


physical decay constant of the radioisotope 
absorption coefficient, cm 
mean standard deviation of cloud, meters 


irradiation time, sec 


gamma flux, shotema/ en ae 


d : - 
ay = horizontal shear, sec 
dv ; = 
cea vertical shear, sec 


Superscript prime (') indicates applicability to cases where radio- 
logical decay is considered 
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APPENDIX B 


DERIVATION OF EQUATIONS 


1. Concentration of a Single Isotope (Instantaneous Release) 


Figure B-1 illustrates a diffusing cloud of fission products 
moving with a mean current velocity v.. At a time t after release, 
the cloud is carried a distance equal to v.t. During this time, 
diffusion and shear effects are causing the cloud to spread and 
elongate. The basic equation for a coordinate system moving with a 


mean current velocity v, is 


2 2 2 
== - (9 y+o2) =, -S4n Sh a O8 : (B1.1) 
y 9x Y oy 9z 


Equation (Bl.1) is a mass balance of the radioisotope concentration 
with neither production nor absorption (losses). The terms on the 
right hand side describe the net leakage of the radioisotopes into a 
unit volume due to diffusion. The second term on the left represents 
net leakage into a unit volume due to a non-uniform mean velocity 
(shear effect). 


Witheanetni tia lconane lon 
S,(O,xs¥»2Z) = M. 6 (x) é(y) 6(z), (Bie2) 


which describes the instantaneous release of the fission product re- 

lease inventory for a single isotope of activity M. from a point source 
ZS : 

placed at x=y=z=0, Carter and Okubo ~ have obtained the following 


solution for infinite space: 
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MeaneeuErent 


flow 


Ditfrusing Cloud 


Initial point Moving with Velocity 


of release 








X 
Y 
Figure B-1. Coordinate System Relationships 
M, 
dl 
S(t, X,y5Z) a a a | ae 
gn 7(a a a tl? 
me) 2 
2 
X 
exp - A 
il 2 2 Zz Z 
ae a+ Fi i «| 
BY 
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ee oy 1 2 2 
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A A 
i ee i ae 
1+ G 2 nt 59 mae 
a x X 22 
> A 2 A, 4 
4A t E + 5 = x + 2, A, met 


t(Q y + QZ )x + = 2 Otay Z 


1 2 v2 a. 2 a 
4A t E vee AT ee 


| 


The concentration of the center of the cloud (x = y = z = 0) is given 


by 


M. 
ae 
i 2am 31/2 A A 
: 1 Z y D> "Zep 
es ) fe \p aa Ts A + QT At 


S (t) = 
Pp 


g7 3/2 


(B1.4) 


If the orientation of the X-axis is in the direction of the mean 
current flow, the concentration as a function of distance along this 


axis is given as 


Mo 2 
S, (t,x) SS ——————— exp 


2 > 
wy ae gt? 4A t (1 + Bt) 


(B1.5) 





2) @©F «4 @@z_=E=p «+ © © Ep & aa &s & 


-—  . s  e ®e —— 
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where 


are (B1.6) 


and 


UL A a 
2 = 1D (Q ry + 2 We (B1.7) 
x x 
Substituting x = x, - vot into Equation (Bl1.5), the following 
equation may be written: 
Mia Ce vet)~ | 
S, (t,x) ee es : (Bi. 8) 


ca \J1 + Bt’ 4A t Glaee Bt”) | 


It is possible to compute from this equation the concentration of a 
Specific isotope at a distance x, measured parallel to the direction 
of mean current flow for an instantaneous release of an isotope of 
activity M, - Radioactive decay has not been considered. 

2. Dose Rate for Instantaneous Release (without decay) 

The dose rate is a function of the gamma flux at the point in 
question. Water has excellent shielding properties and because of the 
high attenuation of the gamma flux, contributions from more than a few 
meters away are negligible. Therefore, it is reasonable to make the 
simplifying assumption that at the point in question there exists an 


essentially infinite cloud of uniform concentration. 
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A person in the water at a distance x,, in the center (x=0) of 
this infinite cloud of uniform concentration S(x,), is exposed to a 


gamma flux as given by Equation (B2.1). 


6G.) = S(xo,t)e "> (A 


Dae 


e + Ave ydv. (B2.1) 


4 


vol. 


Equation (B2.1) comes directly from the fact that the gamma flux from 
an isotropic point source (a volume element dV of concentration S) 


located at a distance r is 


For an infinite cloud, the total flux will be equal to the gamma 


flux due to the volume element integrated over the entire volume. 
Our AUX 
The term (A, e - A, e ) represents the buildup factor. 


Values for the parameters A Ay» Qs: and a, are given by Roelweihias 


Jee 
and are presented in Table B-l. Integrating Equation (B2.1) over 


the solid angle 


Hioxcemt) = S(x,,t) A,e + Ave 


iL i 


—u(1+a, ) xr ~u(1+a,)x 
ais (B2.2) 
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=) Ay 
O(xXo; C) =o ee, ©) —————— naa | (B23) 


u (lta, ) u(1+a,) 


where 6, S, U, Ay» Ay» A> and a, are energy dependent. The dose 
rate in rem per hour for an RBE of 1 due to this concentration is 
equal to the product of the average energy of the gamma photons, 


the gamma flux, and an energy dependent conversion coefficient. 
DICE cers G) or K (E) E Une sc) (B2.4) 


Substituting Equation (B2.3) 


A, (E) A, (E) 
GE. %. 56) =k (CE) ob) S(E.xXe. t) ———_ id 
[u(E) (1ta, (E)) u(E) (1+a,(E)) 
(B225) 


If an energy coefficient, e? is defined by the following equation 


ee A, (E) A, (E) 
e = K(E) EA = = ee ee (B2.6) 
ie u(E) (1+a, (E)) u(E) (lta, (E)) 


|p: 
where ES represents the average energy of the n° energy group of 
a seven group discrete gamma spectrum, then the dose rate due to 


the entire fission product release inventory is 
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/ 
E “t a 
; He a) 


n=1 


Sinan 3 + at? 


(Xo —~ Vo t)* 


DCC.xX5) = exp - (B2.7) 


4A t (1 + Bt’) 


k ai 


3. Total Dose for Instantaneous Release (without decay) 





The determining factor in establishing the exclusion radius 
is the total dose I which is received over a period of time from 
t = 0 to t = T where T is the time required for evacuation of 


personnel from the area. The total dose may be determined from 


T 
I(x.) ra DiGteno) dt, (B3.1) 
O 
i 
d 1 (Xo - vot) 
I(x.) = e M a —— = exp, Soe OE 
a 3/2 2 9 
ee 4: Yjitst 4A tC1+8t) 


(Baq2) 
For determining the low population radius, the applicable 


equation is 


| l 1 (Xo - Hee) 
I(xo) = . e M fo eC (eats 
ee 3/25] 2 2 
! t 1+ Bt 4A t(1+pt ) 
n=1 - x 


(Boo) 





89 


where the upper limit is taken as infinity. 
Equation (B3.2) is best suited for computer solution by 
numerical techniques. An analytical solution is possible for the 


case where the current velocity is zero. For v. = 0, Equation (B3.2) 





becomes 
7 : 2 
1 ae 
xo) = eM 2) eer eee <a dt. 
3/2 see? 4A t(l+st ) 
n=l 
O 
(Boe) 
For t << fe ) (ie =~ 8 hours for typical diffusion coefficients) 
( 7 = : 2 
° 
i.) ok elie ie Bi? expe dt. (B3.5) 
nn | 
=: as 


Z 
Substituting 6 = rs (4A and u = 6/t; du =-(6/t )dt 


we u 3/2 Bena) 
IL Glee ee via € M o' (5) (Guee) _ du 
vgi=sk / - 
6/T 
Ine. a eM? o g 1/2 wl l2 8 du 
n=l 
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This is a form of the incomplete gamma function. 


7 
-1/2 iL | 
LSet a 6 PC > yee 


llosa) = 
| 
= well ~ 
7 
: =1i/2yahee 1 me gplimen 
= eM Sa 6 Co ae oe “at 
n=1 
‘s 
= -1/2 ae 
: eM a 6 \j 7 = yr Erf 7 ; 
=] 
7 
= -1/2 
= ; apes Yu Erte 6/T  , 
nom 
n=l 
or 7 
2a Ae 3) e M 
| x | nn 
ec.) = eee Erfe __ a eee . (B3.6) 
x 2 Ava 
x 
4. Summary - Instantaneous Release without Deca 





The preceding equations have been derived from Carter and 
Okubo's solution in infinite space. No allowance has been made for 
the reflective boundary effect of the ocean floor. As the appli- 
cation for which these equations are intended deals with locations 


at or only slightly above the ocean floor, an approximate solution 





= FE 


may be obtained by multiplying Equation (B1.3) by two, i.e., adding 
the image with respect to the ocean peor The equations below 


include this factor of two. 


Concentration of a specific isotope i; 


S, (5X0 5Vo) ee OND a (B4.1) 
Nise 4A t(1+8t ) 
’ ae -1/2 
Concentration of a specific isotope i at a time t << 8 ; 
sep (xe - vet)” 
S (t5X0sVo) = 2M, at exp - a: (B4.2) 
aan Ss 2 
Concentration of a specific isotope i at a time t >> 8 ; 
Z 
2M,a (Xo — Vot) 
S, (egosvo) = —— oS Oxp —- (B4.3) 
i Sed. 3 
Ct 
/ B 4A Bt 


Dose rate due to the total fission product release inventory; 


ri 
Z e M O 
‘qumal 9 
, (Xo-Vo t) 
n=l 


Dt. ee o.) = Fo) eS (B4.4) 
13/2 “Wise? 4A t(tst ) 
Four hour total direct gamma dose; 
T=4 hr 
Gg sro) = DiGlnee Agee Valte (345) 


4 
O 





O2 


Total dose approximation for zero current velocity; 


iy. 
} e M 
non 
Tse) - SS ____,——- rfc i (Bane) 
eS me 2A T 


55 Decay of Fission Beenaeeee ae 





It is possible to express the rate of decay of the complex 
fission product mixture for the period of time from about 10 seconds 
° -1.2 
to several weeks after shutdown as approximately equal to kt : 
Certainly for times greater than 10° seconds, this approximation is 
conservative as shown by Figure B-2. The following assumptions 


are made: 


(1) The fission product release inventory decays at the same 
rate as the gross fission product inventory. 


(2) Only fission products with a half life greater than 240 
seconds are included in the mixed fission product re- 


lease inventory. 


(3) No decay of fission products occurs until after 1000 
seconds after shutdown. 


(4) After 1000 seconds after SNS CO ny the mixed fission 
product inventory decays as kt ~°* . 


Then the decay of the fission product release inventory is written 
as follows: 
M' = M fomet < ioe seconds, 


i = eh ak ones for ct > 10° seconds. (B5.1) 


6. Instantaneous Release with Decay 


If w(t,xe,v.) is defined by Equation (B6.1) as 





Ie 


10 | ; : 
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a =1.273) 7 \ Qos 6.988 lop 
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aa \ 
gre Ss 
Bs 
Ss 
\ 
=5 l 
10 
10° 10° 10° 10° 10! 
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Figure B-2. Empirical Equations for Rate of Gamma Energy 
Release Following Fission (Data from 
Boderaroc 
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2 
on (Xe=vicits) 
CaeN as Vo) = emai aia aed > 5 © iain 


Z 
| ee 4A t(1+6t ) 


then from Equations (B4.1) and (B5.1) the following equations may be 


(BO. 1) 


written to give the concentration for the total fission product 
release inventory at time t after release: 
| 3 
Sees Ve) = NM (tyke. Ve) for t < 10 seconds (B6.2) 


Meee como” seconds ACbess) 


Brae, Xai, a) 


Similarly, the dose rate may be given as 


DUEexX..Ve) enor & 5 10° seconds, (B6.4) 
-1.2 
& 


it 


Ces Xo ie 


k Dt. xo.ve) pe bomen 10° seconds (B6.5) 


D (ty xoswe@ 


And the total dose received up to time T after release is 


i 
le eee D Guecoe s diteeetor I < 10° seconds (B6.6) 

O 

10° Ar 


alee 


ete ob DCt xO dt ar omic T DiGtrcans oo er 


».. 
~ 


for Ll > 10° seconds. CBG) 
When the concentration of a specific isotope is investigated, 

the appropriate decay constant do must be used rather than Equation 

(B5.1) since Equation (B5.1) is valid only when considering the 


entire mixed fission product inventory. For a specific isotope 


: mee 
S, (t,X0sVe) = Me ‘1'V(t,xesve), ror alle e > Os (B6.8) 
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From Equation (B1.4) it may be seen that the peak concentration 
will be given by 


' M. a e 
ore = (B6.9) 


eis 
when decay is considered. If Sar is set equal to the applicable 
MPCC, Equation (B6.9) may be rearranged as follows: 


2X e Moccmee eZ 
af 


ee Oe Goan). (B6. 10) 


(Bt 


Solution of Equation (B6.10) will give the time at which the peak 


concentration will equal MPCC. 


7. Continuous Release (without decay) 


For the instantaneous release of an amount M2 (where M is the 
fission product release inventory at t=0 and 2 is the leak rate from 


the containment), the concentration along the X-axis will be 





& nee 
Mg [ (Xe Veieo 
Shab og Vo) Sa re 9 9 ee ene dt. (By) 
| oe Vasee, ? 4A to (1+ tc) 
O 
or substituting Equation (B6.1) 
iE 
S.Ct ,. xo oa MUCt> , wes Vo) dG. CBU eZ ) 
4 
O 


For a given distance greater than zero, the concentration will 
be initially equal to zero. If the leak rate is held constant and 


neither depletion of the gross fission product inventory nor 
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radiological decay is considered (both conservative factors), then 
artee a period of time, a steady state value for the concentration 
will be achieved. This steady state value may be achieved by setting 
the upper ine Equation (B7.2) equal to infinity, or 
StXesvo. = Melia. kos Verdes. .« (B7.3) 
O 


The dose rate will be 


co 


Bicomioy = » eM R WiGiers scoop aie . (B7.4) 


n=l 
O 


For steady state conditions, the total dose received will simply 


be the product of the dose rate and the time, or 


7 
MGcov.) = | e M Coie (teks ve) at. (B75) 


44 


n=1 oO 


where T = 4 hours for exclusion radius determinations. 


8. Continuous Release (with decay) 





The equations developed in the preceding section are based on 
two very conservative assumptions: 

(1) steady state conditions are achieved instantaneously, and 

(2) radiological decay is ignored. 

In the case of high velocity currents and relatively short 
distances, actual conditions closely approximate the first assumption. 


Because of this and the fact that such an assumption simplifies 





oe 


computer calculations, it is desirable to retain this assumption and 
base the following equations on the steady state results. 
If decay is considered, the applicable equations at time t may 


be given as follows: 


co 


' 


S (Xo Vo) = | Miu Conxo sve dts. - ome a< 107 seconds (B8.1) 
} 
O 


co 





' — 
S(t 5X0»Vo) = Mgkt noe | ViCEgxes. Ve) deems, tor t > 10> seconds 
Ae 
oo (B8.2) 
' “At f 
S,(t5X0 Vo) = Mte — Wits, Xcgvcl dees, Lor a Specific tcotope 
and all t - 0 O 
(B8.3) 
7 oO 
' ; | 3 
D, (X0»Vo) = eM R | WiGies sees Vo) dite yeh OT ites) 0 
: ret 
seconds, n=1 O (B8.4) 
a co 
: =i, 2 
Dy (t»XosVo) i eM Rkt Ute ea ee el era LOG 
n=1 / 
O 
> 10° seconds, (B8.5) 
ij t 
I, (t,XoVe) =t Dy (Xe sVe)> eo g ee 10° seconds, and (B8.6) 
|e 
r 
g | 7h.2 
I(t, Xo Vo) = (10 sec)D, (Xe5Ve) + D, (Xe sVodk to Giteee as 
bor t - 10° seconds. 107 


Ceca } 





98 


The steady state concentration as a function of distance is shown in 


Figure B-3. 


Log Concentration, S&S 






S(x.), no decay 






S'(x.), after decay time 


th >10~ seconds 






) NX 
SS € oe nl ion XN 


fime 6, > et 
Z 


Log Distance, X, 


Figure B-3. Concentration of Fission Products 


as a Function of Distance 


If the instantaneous steady state assumption is not used, the 


results may be summarized as follows for times after release t greater 


than 10° seconds: 


S(t sX0,Vo) = Mekt_ 


and 


1s 


a2 
. | PiCtiog emive) Gee 


O 


i 
: a: e 
Dy (tsXo»Vo) = / eM Rkt Wilton Xo sevice 


(B8.8) 


(B8.9) 
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' ' 
where 35 and D, are the concentration and dose rate respectively 


at time t with radiological decay considered. The total dose received 


for exposure times in excess of 1000 seconds is 


10° t 
j t t 
I, (t»Xo Vo) = D, (Xo »Vo)dto + D,(tosXo»Vo)dto (B8.10) 
O 103 
5 3 
; 10 10 
I, (tsXo»Vo) = | eM, Q Gites ||. 4) (tae cere eed te, 
n=l O O 
ie ie 
-1.2 
+ Kea de, |WGhs. X05 Ven Glesel. (B8.11) 
10° 0 


9. Continuous Release - Approximate Analytical Solution for v,.=0 


t 
Factoring D, (Xo »Vo) from Equation (B8.7) and substituting Equation 


(B8.4), it is possible to obtain, for t > 10> seconds, 


ree t 

7 | 

' 3 -1.2 
I(t »Xo»Vo) = eM Q itt, 10°+k] t, dt 
n=l 10° 
(B9.1) 
ae 7 -1/2 

Dividing the integral over  W(t.,xX.,v.)dt, at t.=B and 


rearranging terms 
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t = 
| | 
I,(t,Xo0,Vo) = 10° + k ee Gite eM Q 


j 
10° nel 


ES = | 


(Co sxoswe dts = ) (te semen cies ° (B9.2) 


Substituting Equation (B6.1) for ~ and making the approximation 
that for times less than t, the term Gece. ol and that for times 


greater than ts the term (148e2)e ete it is possible to obtain 


f ~ 


i 

t = 
L2GChj%e.V5) = 10° +k to * dt. x M 220 

2 nm, 

10° n=1 _ 
Ee on 
= 2 2 
-3/2 (Xo-Vot.) sie, =5/0 CRASS), 
xen 4A Bt. 
xX 
Lo c 
= i 
(B9.3) 


For the case where v,.=0 and proceeding in a manner similar to 


Section 3, the following results are obtained: 





om 


/ 


t 
I, (t,Xo) = 10° + 5k(10 Ue 220, ys M 
i 71 


2 yma, 
= - PEPE 
Xo 
{ 3 = 
I, (t,x0) = |10° + 5k(10 
Xo 


ErEe 


24/4 Ge 
a: 


HO. Virtual Distance 


+ 


uF 


7 
42a 4/7A 
6 : ene 2, x \s M 
Xe nn 
n=l 
i ac 
3 Erf —_—__—_ : (B9.4) 
2A: 
x i 


A virtual distance concept is included in the computer solutions 


in order to examine distances very close to the containment. If the 


effective area of the diffusing cloud is defined as extending a 


e g Z o & 
distance of 30, where for small times, o is given by Carter and 


olen yee as 


2 See E1051) 


pan mete . (B10. 2) 
xy 
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then the cross-sectional area of the cloud (Y-Z plane) may be approxi- 
mated as 


A = 91 cor , (BO. 3) 


IS, io 


A’ : B10.4 
; ( ) 


> 
tt 


Let ce A 
x 


If the area of the cloud is set equal to the area of the containment 
the time required for the diffusion process to approximate actual 


conditions will be given by 


Sa (B10.5) 


or for a given current velocity, the virtual distance will be 


Z 


Wigecl 


N ee : (BIG) 
72 At? ae’? 
x 7 


Figure (B-4) illustrates the virtual distance concept. 


m1, Investigation of the Dose Level in the Vertical Direction 


In order to examine the dose level as a function of the 
distance above the point of release for an instantaneous release 
case without decay, it is necessary to begin with Equation (Bl.3). 
If x and y are set equal to zero and it is assumed that the mean 
current velocity in the z-direction is zero, then this equation 


becomes 
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A A 
lige 2 ey), Wek? Wize 
=3/2 tts ty a t3% wt, 
oy lal ie x x 
SC 9 Maes meee (00) een 5 
(B11.1) 


Allowing for the reflection from the sea floor, the dose rate 


due to the mixed fission product release inventory will be 


7 
2 as A A 
1 + (Lp? 24194 2 )-? 
— 12 yA Bes A 
n=l x x 2 
i) ek — 5 Zz 
3/2 1/4 2 gt? 4A t( 1+ et) 
(B11.2) 
current 
= 
velocity, vo 
reactor cc: 
containment ec ee 
SS Sy ' ‘ ‘ 
JI 1 > 
virtual distance xX, 
distance | 


Figure B-4. Virtual Distance Concept 
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APPENDIX C 


DETERMINATION OF THE MAXIMUM PERMISSIBLE CONCENTRATION 
OF A RADIOISOTOPE IN SEAWATER (MPcc)14 


Case I. See rene VAG Lyk eEMod 


For those cases in which the limiting intake is determined by 


body burden, then 


Sa 
Meee = ee (Ce yy (C1.1) 
nf B 
nb 


MPCC's for radioisotopes in Case I are given in Table C-l. 


Case II. Concentration Factor Method 
For those cases in which the intake is limited by irradiation 


of the gastrointestinal tract by its contents, then 


SMe G) ly 
ME CC = tar ae ee CCl) 


it 
ns 


where F is the ratio of the consumption of water, in liters per 


kilogram, to the consumption of seafood. The factor F is approxi- 
I 
S ‘ : : 
mately equal to 10. The term — is equivalent to the concentration 
nw 


factor of the seafood, f. Hence, Equation (C1.2) may be written 


S (MPC) F 


F (Cie3) 


MPCC 


MPCC's for radioisotopes in Case II are given in Table C-2. 
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The natural chemical abundance of the carrier in seawater is an 
extremely important factor in calculating MPCC's. This factor enters 
into the calculations either directly (specific activity method) or 
indirectly through the concentration factor (concentration factor 
method). Table C-3 gives the average chemical abundances of the 
elements in seawater. ol ontae points out that such a tabulation as 
this suffers from the fact that many of the values are from a single 
set of analyses of surface waters and as such are quite possibly 


uncharacteristic of the ocean environment as a whole. 
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Calculations of MPCC for Radioisotopes for Which Limiting 


Intake is Determined by Irradiation of the Gastrointestinal 


Concentration Factor Method 


Critical 
Organ 


GI(LLI) 
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GI (ULI) 
GUCEED) 
GI (ULI) 
GI (LLI) 
GI (LLI) 
GI (LLI) 
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Table C-2 (continued) 
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MPCC 


f MPC 
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APPENDIX D 


FISSION PRODUCT INVENTORY 


The fission product inventory to which the release assumptions 
of Chapter 4 are applied in order to obtain the fission product release 
inventory is presented in Table D-l. Table D-1l includes all fission 
product radioisotopes with a half life of over four minutes (6.67 x noe 
hours) which decay by the emission of gamma photons of over 0.1 Mev 
energy. In addition, the inventory includes those fission products 
which are not gamma emitters but for which maximum permissible body 
burdens have been established. 

The inventories are calculated, using ORNL Dy on the basis 
of an irradiation time of 3 x 10’ seconds (approximately 347 days); 
time after shutdown equal to zero; and a neutron flux of en neutrons 
per square centimeter per second. 


Solubility informaton was obtained from the Handbook of Chemistry 
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Table D-l. Fission Product Inventory 


Element Isotope Fission Product Inventory Half Life Semup a ti ty 
in 
By Element By Isotope Cold Water 
(curies /MW (curies/MW) (hours) (mg/1) 
Zine < 1 
| 1 ; 
72 a e290 x 10 insol. 
Gallium it insol. 
; an 
72 —. L 1.43 x 10, 
73 ak 500° a0 


Germanium 191 insol. 
75 7 1.37 x 10° 
a7] 52 eZ Os x 10, 
78 2 ao xe 
Arsenic IOUS Onis Onl 
1a) i. Sy oie aes 10° 
78 ley 2 oie | x 101 
79 300 iO. x 104 
81 170 67 xo 
Selenium dog Lore insol. 
79 <1 5.69 x 10°, 
81 75 Zoe x 104 
81m 5 OFA? x 104 
83 ie 1) fle Eelp 
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Table D-1. (continued) 


Element Isotope Fission Product Inventory Half Life Solubility 
ala 
By Element By Isotope Cold Water 
(curies/MW) (curies/MW) (hours) (mg/1) 
Bromine 1 320 6%50 x 107 
82 10 3.59 x 10° 
83 4,050 Ze40 x 10-5 
84 9,460 OU x 0 
Krypton 712.209 IL SOD 
83m 4,050 1.90 x 10° 
85 119 9.00 x 10, 
85m [22500 deeb NS) x 10 
87 259300 ole SHOP 5s 10) 
88 32,200 Zeal? x=lo 
Rubidium Zee 26 decomp. 
86 26 4.68 x 10° 
88 31,900 2396 & 105 
89 39200 Vacs) Wie 105 
91 505500 2204 xO 
Strontium 795.3530 decomp. 
89 37,200 1.30 x 10; 
90 T2540) ee SX 105 
ell 50e2800 oi & 10) 
92 53,000 Zan x 10.5 
93 54,100 Laie x eo 
Metrum 291,740 sl. decomp. 
90 1,340 6.45 x 10, 
91 51, 000 las x 10°75 
91m 20,200 oh 10s 10, 
92 S25 (410,0) B00 x 10, 
93 55.000 W000! x 10.5 
94 565200 As Same 105 
95 55,600 EwS 2-10 
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Table D-l. (continued) 


Element Isotope Fission Product Inventory Half Life Solubility 


ei 
“By Element By Isotope Cold Water 
(curies/MW) (curies/MW (hours) (mg/1) 
Zirconium 1067300 3 insol. 
95 525200 loa 10, 
97 53,800 LO x ako 
Niobium JA 24 i insol. 
93m 72500 Bale x 10, 
95 51,100 Oea40 "x 10, 
95m 519 9.00 x 10) 
96 5 2.58 x 10, 
sii 52,400 1276 x v0 
Molybdenum VZ6—7 00 l insol. 
99 51,400 6.2/0 xX 10°71 
ed fk SHUG V1 (EE Se 105 
102 36,000 22.00 x a0 
Technetium 46,780 0 ? 
99m Dy low 6.04 x 10°43 
TO) 1010 oS Xone 
Ruthenium 257540 9 insol. 
103 24,600 O.ca x 10, 
105 7,500 2A 0 x 10. 
106 1,520 8.76 10°, 
107 1,720 6567 xO 
Rhodium OSe 722 # insol. 
LOS PG E1020, 0700 x 10, 
105 7-560 OFiO Dx 10.5 
£07 ge 770, Ae BY ore 10, 
109 242 P2000 x 210 
Lead 502 1 insol. 
109 255 1236 = 10" 
jill 16 &. 67 x 10, 


ee? oA. Z2e80 x LO 











Element 


Silver 


Cadmium 


Indium 


Tin 


Antimony 


Isotope 


110m 
111 
gL 
143 
HES: 


iil s: 
HLS: 
JL i 
AWM faye 
118 


115m 
7 
ee 
118 
119 


117m 
119m 
al 
ZS 
TZ 
125 
126 
1 


i225 
126 
127 
128 
ae 
130 
aie iik 


Table D-l. (continued) 


Fission Product Inventory 


By Element By Isotope 


(curies/MW) (curies/MW) 


460 
30 
Loy 
95 
8/7 
91 


350 
87 


86 


87 


366 


35202 
come hl 
eal 
iLike: 


hd eee) 
144 
862 
2,080 


55,597 
41 
856 
DAW) 
4,300 
8,860 
17,200 
22,200 


Is 


Half Life Solubility 


aie 
Cold Water 
(hours) (mg/1) 
3 insol. 
6.48 x 10, 
eo2. xX 10, 
320) xX 10, 
aeaO x 10-1 
3233. x 10 
l insol. 
550) x 10, 
IOs: x 10°13 
Sag5 x 10, 
Vd SHO) ape 1075 
5007 x 0 
0 insol. 
Ao) x 10, 
JeatO: x 10, 
1.90 x 10, 
7.50 x 105 
2.92 x 10 
D ab oh syepke. 
3a30° xX 10. 
6,60 x 10, 
234 3: xX 10, 
Oeer x 10°11 
6.59 x 10, 
Zo 105 
Sas. x 10, 
iO x eke 
4 insol. 
Zo Xx 10, 
9.00 x 10, 
GEO x 10, 
lO) x 10, 
Ae60 x 1075 
200) x 104 
Seaun x 2L0 


stege 








7.4.2. 


dal 


Table D-l. (continued) 


Element Isotope Fission Product Inventory Half Life Solubility 


ii} 
By Element By Isotope Cold Water 
(curies/MW) (curies/MW) (hours) (mg/1) 
Tellurium 187,143 3 insol. 
125m 5; 1.39 x 10) 
12 2,030 9.30 x 10, 
127m 278 2lGy x 10) 
129 8,350 12 Ox 10, 
129m 2,890 ioe xX 10°54 
Isa 24,400 4,14 x 10, 
131m 3,890 3,00) x 10; 
182 373300 7.78 x 10, 
133m 51, 800 Tee: x 10/1 
134 56,200 ime x 10 
Iodine 233,800 9 2597 us 107 
31 23,600 ie93) x 10, 
ie 2 3535100 20 x 10, 
133 55,100 2708 xX 10°, 
134 65,400 Creo 1X 10, 
NSS 51,600 6.68 x 10 
Xenon 125,368 » 241 cm /e 
131m 258 2.08) x 10, 
iss 1,350 1.26 x 10) 
Son 56,200 5 oe x 10) 
135 3,760 O20 ox 1075 
135m 15,400 2.00. x 105 
138 48,400 2465 1X LO 
Cesium 218) Ja 9 decomp. 
136 212 OL Zax 10, 
137 L026 Zao x 10°75 
138 Ne 5018 Diets 0.x 105 
i? 49,200 Uasyoimeg U8, 
Barium 203,000 0 decomp. 
139 52,400 1.42 x 10, 
140 54,000 Sm0 7 x 1074 
LA 49,500 S200: x 103 
142 47.100 1.00 x 10 
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Table D-l. (continued) 


Element Isotope Fission Product Inventory Half Life Solubility 


in 
By Element By Isotope Cold Water 
(curies/MW) (curies/MW) (hours) (mg/1) 
Lanthanum 1 decomp. 
140 54 ,000 4.02 x 10, 
141 51,000 3 Oa 10, 
142 51,000 LZ x 1071 
L428 54,100 3.17 x 10 
Cerium 210,100 D insol. 
141 51,800 Tos. x 10, 
143 53,500 a0 x 10, 
144 29,200 6.96 x 10°71 
146 27,100 2032 x LO 
Praseodymium 146,000 D decomp. 
143 a2 00 3.29 x 10/1 
144 28,700 2 oD x 10, 
145 36,200 6.00 x 1075 
146 29,000 4.06 x 10 
Neodymium is ;soL 9 decomp. 
147 Zoe Za x 10, 
149 11,500 2 U0nx 1074 
151 4,140 Jarra Oe ALO 
Promethium 19,970 h ? 
147 4,210 ISI. oie 10, 
149 11,600 5.40 x 10) 
151 4,160 Zit OUR LO 
Samarium Zao 5 insol. 
isl ome SS x 10, 
153 Paps) 4.70 x 10°13 
15> 271 Bien x 10, 
156 T2 1.00 x 10 





118 


Table D-l. (continued) 


Element Isotope Fission Product Inventory Half Life Sodas lay 


in 
By Element By Isotope Cold Water 
(curies/MW) (curies/MW) (hours) (mg /1) 
Europium 455 r insol. 
155 3 1 RS Jo aa 10, 
156 SZ. 3/7 0ux 10, 
JLo, 63 ee x 10, 
lisyss iy, 1.00 x 10 
Gadolinium 10 1l insol. 
159 10 Laon 10 
Terbium <ome ” insol. 


161 al 1.68 x 10 


Mh 








J 
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APPENDIX E 


"SEADIF' COMPUTER PROGRAM 


Two separate and complete computer programs have been written. 


SEADIF I (Table E-1) is a FORTRAN IV computer program for the uncon- 


tained reactor system whereas SEADIF II (Table E-2) treats the con- 


tained system. With the information shown in Table E-3, provided in 


the format given by Figure E-l, the following results may be obtained: 


a. 


Dose rate as a function of the distance from the reactor 
at four hours after release 

The integrated dose received in the first four hours after 
release as a function of distance from the reactor 

The integrated dose received for an infinite period 

of time after release as a function of distance from the 
reactor 

The exclusion radius 

The low population radius 

The time (or in the case of SEADIF II - the distance) 
required for the concentration of a specific fission 


product to reach MPCC. 


SEADIF I evaluates Equations (B6.4), (B6.5), (B6.6), and 


(B6.7) in determining the dose rate and total integrated dose 


received. 


In the low population dose, an upper limit is set by 


successively evaluating Equation (B6.5) until the numerical value of 


EhiSmitiwetven is lessmthan lWa@00 themvalwme o@& thesiwaction at its 


maximum. 
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SEADIF II evaluates Equations (B8.1) through (B8.7) by numerical 
integration techniques employing the Gaussian numerical integration 
library function DGNTGR. 

Both programs require that current velocity and distance from 
point of release be provided as input data. The programs provide 
complete exclusion radius, low population radius, and dose rate 
results for each current velocity (up to 16 maximum on any one run) 
and safety radius results for the maximum current velocity. A zero 
velocity input may be used although such an input will require more 
computer time, especially for the continuous release case. 

It is important that the distance used be ordered smallest to 
largest and the interval between distance points be kept reasonable. 
Reasonable distance intervals (on the order of 500 meters) result in 
computer time savings. Under no conditions should the distance 
intervals exceed 5,000 meters. Up to 40 distances from point of 
release may be provided as input data. The minimum distance that may 
be used is five meters. However, because of fundamental release and 
decay assumptions, results for distances of less than 20 meters must 
be used judiciously. At these small distances, round off error from 
the numerical integration techniques may become large. Excessive 
roundoff error will be noted by a numeral 1 appearing under the 
column headed ERROR in the program output. 

Both programs have a virtual distance option. Providing the 
containment diameter on the appropriate card in Table E-3 will correct 
all distances by an amount equal to the virtual distance, V, as 


computed from Equation (10.6). A zero containment diameter or a blank 


“<—p 
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card will call for a zero correction. Cards for up to 200 isotopes 
may be provided as input for the fission product release inventory 
calculations. For the results presented in Chapter 6, the data shown 
in Appendices C and D was used. Omission of any isotope results in 
a smaller release inventory and corresponding reductions in the dose 
rates and distance factors. It is important to include all fission 
products with a half life of over four minutes and their decay 
information by energy group whether or not an MPCC has been 
established. Omission of the MPCC from a data card results in the 
omission of the peak concentration (safety radius) calculations for 
that isotope and does not effect the release inventory. 

If it is desired to change the release percentages for the 
soluble and insoluble isotopes, or for that matter, any single 
isotope, this may be done in the last data field of the isotope 
information card. 

Representative computer time for five velocity cases and 15 
distances is about 30 seconds for SEADIF I and about 100 seconds 
for SEADIF II. For SEADIF II, the smaller current velocities require 


more computer time. 
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Table E-l. FORTRAN IV Computer Program for an 
Uncontained Reactor System 


KAKKKSEADTP**** 


INSTANTANEOUS RELEASE CASE 
WITH DECAY 


IMPLICIT REAL*8(A-H,0-Z) 
EXTERNAL FCT,WOD 
COMMON ALPHA,BETA,AX ,NX/ISOT/IDENT(200) ,AMT(200) ,PCC(200), 
IDIST(200) ,TIME(200) ,PER(200) ,HALF(200)/FUN/VEL, XX 
DIMENSION VO(16) ,X0(40) ,GAM(200,7),A(200,7) ,ANS(10) ,POP(40) , 
1SUM(7) ,EN(7) ,TEMP(7) ,AUX(1000) ,RINT(40) ,DOSER(40) , VDIS(16) 
PRINT 10 
10 FORMAT ('1*****SEADIF***** INSTANTANEOUS RELEASE CASE‘) 
READ 11,POWER 
11 FORMAT (F10.4) 
PRINT 12,POWER 
12 FORMAT ('OPOWER=',F10.4,' MW') 


READ IN CURRENT INFO AND COMPUTE ALPHA AND BETA 


READ 13,AX,AY ,AZ,OMEGAY ,OMEGAZ 
13 FORMAT (2F10.3,3D10.3) 
ALPHA=1 .0/8.0/3.1415927**1 .5/DSORT (AX*AY*AZ) 
BETA=( OME GAY **2*AY/ AX+OMEGAZ**2* AZ/AX)/12.0 
PRINT 14,AX,AY,AZ,OMEGAY ,OMEGAZ,ALPHA,BETA 
14 FORMAT(' AX=',F10.5' METERS SQUARED PER SEC'/ 
leave El0e5,”~—C METERS#SQUARED PER SEC'/ 
2' AZ=',F10.5,' METERS SQUARED PER SEC'/ 
3 BUmeGh Y=" 010.3,° PER. SEC”/ 
4" OMEGA z=" ,DI0.3, PER ieee, 
5' ALPHA =',D10. 3, SEC*3/ GaEeER METER*3 / 
6° BETA =' ,D10. 3,' PER SEC SQUARED') 
READ 15, NV, NX 
15 FORMAT (213) 
READ 16, Mae joa 
READ 17,(X0(I),1 
16 FORMAT (8F10. 4) 
17 FORMAT (8F10.2) 


= Le 
=1,NX) 


> 


READ IN EQUIVALENT CONTAINMENT DIAMETER (METERS) AND COMPUTE THE 
VIRTUAL DISTANCES FOR THE GIVEN VELOCITIES. 


READ 18,DIAM 
18 FORMAT (F10.2) 
PRINT 19,DIAM 
19 FORMAT (' CONTAINMENT DIAMETER =',F10.2,' METERS') 
VIRT=DIAM**2/7 2.0/AX**0.75/AY**0. 25 
DO 20 I-],NV 
20 VDIS(1)=1)*VIRT 


4 — 
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Table E-l. (continued) 


READ IN ISOTOPE INFO AND COMPUTE SOURCE TERM 


READ 21,NI 
21 FORMAT () 3) 
READ 22,((IDENT(I) ,AMT(1),PCC(1I) ,HALF(1),(GAM(I,J),J=1,7), 
IPER(I)),1=1,NI) 
CONFORM EM (15 aes O23 7F5. 1,65. 3) 
DO: 22> 1=|.,NI 
De 23-0217 
A(1,J)=AMT (1)*GAM(1,J)*PER(1)*0.01*POWER 
23 CONTINUE 
DO 24 J=l1,/7 
SUM(J)=0.0 
DO 24 [=1],NI 
SUM(J)=SUM(J)+A(T,J) 
24 CONTINUE 
EN(1)=4.30D-15 
EN(2)=1.46D-14 
EN(3)=2.04D-14 
EN(4)=3.17D-14 
EN(5)=2.05D-14 
EN(6)=4.17D-14 
EN(7)=5.26D-14 
XEM=0.0 
DO 25 J=l,/ 


TEMP(J)=EN(J)*SUM(J) 
XEM=XEM+TEMP( J) 
25 CONTINUE 
XEM=XEM*3./7D+10 
PRINT 26,XEM 
26 FORMAT (' XEM='D10.3,' RM*3/SEC') 


DETERMINE DOSE RATE AT 4 HRS AND TOTAL DOSE RECEIVED IN 
FIRST 4 HRS. 


PRINT 28 
28 FORMAT (']IRADIATION DOSAGE INFORMATION'/10X,'DOSE RATE IN' 

1' R/HR AT TIME=4HRS AFTER RELEASE'/10X,'GAMMA DOSE ' 
2'RECEIVED UP TO 4 HRS AFTER RELEASE'/10X, ‘EXCLUSION RADIUS' 
3' BASED ON DOSE NOT TO EXCEED 25R'//'O DIST(M) VEL(M)' 
4'SEC) DOSE RATE(R/HR) DOSE(R) LOW POP DOSE 
5" ERROR LITO INTERVAL VIRT DIST(M)') 

EPS=0 .01 

NDIM=12 

DO 64 I[=1,NV 

VEL=VO(1) 

DO 41 J=1,NX 

TRY1]=0.0 

TRY2=0.0 

XX=X0(J)+VDIS(1) 


? 
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3 
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32 


33 


34 


So 
36 
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38 


32 
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Table E-l1, (continued) 


DRT=1 .44D+4 

DOSER(J)=XEM*ALPHA*FCT( DRT) *2.86632D+7 

CALL LIMIT(XX,VEL,XL,XU,IERR) 

IF (XU. GT.1.44D+04)XU=1 .44D+04 
(VEL.LT.0.02) GO TO 29 

F (XU.EQ.1.440+04 ) XL=L=1 .00D+03 
(XL.EQ.1.00D+03) IERR=5 

[pe RET 1000.0) “GOrTO0 37 

IF (XL.GE.1.44D+04 ) XL=1.00D+03 

CALL DQATR(XL,XU,EPS,NDIM,FCT,S, IER, AUX) 

RINT(J)=XEM*ALPHA*S*7.962D+3 


DETERMINE LOW POPULATION DOSE 


SAVE=0 .0 

K=1] 

XL=XU 

XU=2 .0*XU 

IF (XL.LT.1000.0) GO TO 39 

CALL DQATR(XL,XU,EPS,NDIM,FCT,ANS(K) ,IR,AUX) 
PEXTR=XEM*AL PHA*ANS (K)*7 .962D+3 
SAVE=SAVE+PEXTR 

IF (RINT(J).EQ.0.0) GO TO 34 

IF (PEXTR/RINT(J).LE.0.01) GO TO 35 

IF {(K.E0.10) GO.10 35 

K=K+] 

GO T0 32 
POP(J)=RINT (J) 
PRINT 36 »X0(J) 
FORMAT (! ar) 
GO TO 4] 
XIT=1000.0 
TRY2=0.0 
IERR=6 

CALL LOWER(XX,VEL,XLW,XUW, IEW) 

CALL DOQATR(XLW,XOW,EPS,NDIM,WOD,SADD, IER, AUX) 
TRY 1=XFM*ALPHA*SADD*2 .0 

IF (XUWW.LT.999.9) GO TO 38 

[PE OCGRO) GO 10) 32 

CALL DQATR(XI,XU,EPS,NDIM, FCT,S,IER2,AUX) 
Pe ,uveeEPiA Ss” / .902D43 

RINT(2)=TRY1+TRY2 

TER=IER+IER2 

GO T0 31 

XI=1000.0 

IF (XU.LT.1000.0)XI=XU 

TERR=IERR+10 

PALSO=0.0 

CALL DQATR(XL,X1,EPS,NDIM,WOD,ANS(K),1R,AUX) 
PADD=XEM*AL PHA*ANS (K) *2.0 

IF (4U.LE.1000.0) GO TQ 40 


+SAVE 
, VEL ,DOSER(J) ,RINT(J } »F POP(J), IER, IERR,K,VDIS(I) 
Decne PO 4 SxS. 4.5%) .2( 12,8 12, 5X,F10.5) 
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Table E-l. (continued) 


CALL DOQATR(XI,XU,EPS,NDIM,FCT,ADD,IR,AUX) 
PALSO=XEM*ALPHA*ADD*7 .962D+3 

40 PEXTR=PADD+PALSO 
GO TO 33 

4] CONTINUE 


DETERMINE EXCLUSION RADIUS BY EVALUATING INTEGRAL AT IOM 
INTERVALS UNTIL VALUE OF INTEGRAL IS LESS THAN 25R. 


IF (RINT(1)-25.0) 42,42,44 
42 PRINT 43,V0(1),X0(1) 
AS FOR Me awk VEROCIIVAOF 3610.4." sM/ SEG. THE EXCLUSIONS 
1' RADIUS IS LESS THAN',F10.3,' METERS') 
GO TO 56 
44 DO 54 K=2,NX 
IF (RINT(K)-25.0) 45,54,54 
45 XX=X0(K-1)+10.0+VDIS(1) 
DO 52 L=1,500 
CALL LIMIT(XX,VEL,XL,XU,IERR) 
IF (XU.GT.1.44D+04 ) XU=1 .44D+04 
IF (VEL.LT.0.02) GO TO 46 
IF (XU.EQ.1.44D+04 ) XL=1.00D+03 
46 IF (XL.LT.1000.0) GO TO 50 
CALL DQATR(XL,XU,EPS,NDIM,FCT,S,IER,AUX) 
XCLU=XEM*ALPHA*S*7 .962D+3 
47 IF (XCLU-25.0) 48,49,49 
48 XXX=XX-VDIS(1) 
PRINT 43,VO(I) ,XXx 
G0) Oh leis 
49 XX=XX+10.0 
GOmnO 52 
50 XI=1000.0 
TRY2=0.0 
CALL LOWER(XX,VEL,XLW, XUW, IEW) 
CALL DOATR(XLW,XUW,EPS,NDIM,WOD,S ,IER,AUX) 
TRY 1=XEM*AL PHA*S*2 . 0 
IF (XUW.LT.999.9) GO TO 5] 
IF (XU.LE.1000.0) GO TO 51 
CALL DQATR(XI,XU,EPS ,NDIM, FCT ,SADD, IER2, AUX) 
TRY 2=XEM*ALPHA*SADD*7 .962D+3 
51 XCLU=TRYI+TRY2 
GO TO 47 
52 CONTINUE 
PRINT 53,VO(I) 
53 FORMAT ('OFOR VELOCITY OF ',F10.4,' M/SEC, AN ERROR EXISTS' 
1'REGARDING EXCLUSION RADIUS') 
54 CONTINUE 
PRINT 55,VO(I),X0(NX) 
SS enenoe OROR VELOCITY OF *,FIl0.4,. M/SEC. THE EXCLUSION ~ 
1'RADIUS IS GREATER THAN',F10.3,' METERS') 
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Table E-1 (continued) 


FIND LOW POPULATION RADIUS(APPROXIMATE METHOD) 


56 TpeePoeen) 25.0) 57,57,59 
57 PRINT 58,X0(1) 
58 FORMAT (' ',35X,'THE LOW POPULATION RADIUS IS LESS THAN',F10.3,' 
IMETERS'/) 
GO TO 64 
99 DO 60 N=2,NX 
IF (POP(N).LT.25.0) GO TO 62 
60 CONTINUE 
PRINT 61,X0(NX) 
61 FORMAT (' ',35X,'THE LOW POPULATION RADIUS IS GREATER THAN' 


135! CRe METERS 7”) 


GO TO 64 

62 XCRE=(POP(N-1)-POP(N)/(XO(N)-XO(N-1)) 
PLUS=(POP(N-1)-25.0)/XCRE 
RAD=XO(N-1)+PLUS 


PRINT 63,RAD 

63 FORMAT (' ',35X,'THE LOW POPULATION RADIUS IS APPROXIMATELY ' 
| splleeese MEPERS 7) 

64 CONTINUE 


PEAK CONCENTRATION LEVELS 


PRINT 66 
66 FORMAT ('ICRITICAL ISOTOPE INFORMATION'//' ID'.8X,'CURIES/Mw' 
1,5X,'MPCC(UC/KG TIME (HRS) DISTANCE(M) HALF' 
2'LIFE(HRS)") 
VV=VO(NV) 
CALL PEAK(VV,NI,ALPHA,BETA, POWER) 
STOP 
END 


SUBROUTINE PEAK(VV NI ,ALPHA,BETA,POWER) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON /ISOT/IDENT(200) ,AMT(200) ,PCC(200) ,DIST(200) , TIME(200) , 
1PER(200) ,HALF(200) 

BE TAH=BETA*1 2 .960+6 

DO 113 I=1,NI 

Mgect 1) .LE.0.0)TIME(1)=0.0 

ie (nee EE.0.0) GO 70 113 
COF=(AMT(I)*ALPHA*POWER*PER(I)/PCC(1)**2*8.573388D-5 
XLAM=0.693/HALF (1) 

KOUN T=0 

T=1000.0 
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Table E-l1. (continued) 


100 IF (BETAH*T**5+T**3.GE.COF*DEXP(-2-..0*XLAM*T)) GO TO 101 
IF (KOUNT+1.GE.100) GO TO 109 
T=T+1000.0 
GO TO 100 
101 T=T-900.0 
KOUNT=0 
102 IF (BETAH*T**5+T**3.GE.COF*DEXP(-2.0*XLAM*T)) GO TO 103 
IF (KOUNT+1.GE.10) GO TO 109 
T=T+100.0 
GO TO 102 
103 Tele ene.. 
KOUNT=0 
104 IF (BETAH*1**5+1T**3.GE .COF*DEXP(=2.0*XLAM*T )) GO TO 105 
IF (KOUNT+1.GE.10) GO TO 109 
T=T+10.0 
GO TO 104 
105 T=T-9.0 
KOUNT=0 
16 Ween *onl~*35Ge-C0r VEXP(-2.0*XLAM*T)) GO TO 107 
IF (KOUNT+1.GE10) GO TO 109 
=aeells, O 
GO TO 106 
107 TE GiGi 1GG00 G0 m0. 111 
T=T-0.9 
KOUNT=0 
108 IF (BETAH*T**5+T**3.GE.COF*DEXP(-2.0*XLAM*T)) GO TO 111 
IF (KOUNT+1.GE.10) GO TO 109 
T=T+0.1 
GO TO 108 
109 PRINT 110, IDENT(I) 
110 FORMAT (' ',15,7X,'**ERROR MESSAGE-T GREATER THAN 1.0D+05 ' 


IRS) 

TIME(1)=0.0 

GO TO 113 
111 TIME(I)=T 

DIST(I)=TIME(1)*VV*3.6D+3 

PRINT 112, I1DENT(I) ,AMT(I) ,PCC(I),TIME(I) ,DIST(I) ,HALF(TI) 
12 5RoRvam ¢* *,15,5(5x5010.. 3) ) 


113 CONTINUE 
CRIT=TIME(1) 
IC=IDENT(1) 
DO 114 [=2,NI 
IF (TIME(1I).GT.CRIT)CRIT=TIME(1) 
IF (CRIT.EQ.TIME(1I))C=IDENT(I) 
114 CONTINUE 
PRON 115, TC VV 
115 FORMAT ('OTHE CRITICAL ISOTOPE IS ',15,'. DISTANCES FOR ' 
WVELOCINY=",F10.4,' M/SEC’) 
RETURN 
END 
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Table E-1. (continued) 


SUBROUTINE LIMIT(X,V,XL,XU,IERR) 

LIMIT ERROR MESSAGE HAS THE FOLLOWING MEANINGS: 

-~O- UPPER AND LOWER LIMIT LESS THAN 4 HRS. NO ERROR 

~]- UNABLE TO ACHIEVE DESIRED ACCURACY ON UPPER LIMIT. 

-2- UNABLE TO ACHIEVE DESIRED ACCURACY ON LOWER LIMIT 

-3- UNABLE TO ACHIEVE DESIRED ACCURACY ON EITHER UPPER 
OR LOWER LIMITS. 

-4- UPPER LIMIT SET EQUAL TO 4 HRS 

~5- UPPER LIMIT SET EQUAL TO 4 HRS, LOWER LIMIT SET EQUAL 
TO 1000 SECONDS. 

IMPLICIT REAL*8(A-H,0-Z) 

IF (V.EQ.0.0) GO TO 206 

IF (V.LT.0.01) GO TO 205 

IERR=0 

1=0 

TSET=X/V 

T=TSET 

TI=0.5*TSET 

IF (V.GT.0.149)TI=0.25*TSET 

TMAX=FCT(T) 


Valea 


Il=I+T 

IF (I.£Q.20) GO TO 204 

IF (FCT(T)/TMAX.GE.0.001) GO TO 200 
XU=T 

I=0 

TL=TSET 

IF (V.GT.0.149) GO TO 203 
Me= ie 2Z..0 

I=I+] 

IF (I.£Q.20) GO TO 207 

IF (FCT(TL)/TMAX.GE.0.001) GO TO 202 
XL=TL 

GO TO 208 

e=On75 IL 

[=I+] 

IF (I.EQ.20) GO TO 207 

IF (FCT(TL)/TMAX.GE.0.001) GO TO 203 
XL=TL 

GO TO 208 

XO=11.0*TSET 

IERR=1 

GO TO 201 

XU=1.44D+4 

XL=X 

IF (XL.GT.XU)XL=XU 

IERR=4 

GO TO 208 

XU=X*] .OD+6 

XL=X 

GO TO 208 

KE=1 07 2-7 0**20*1SET 





208 


300 


301 


302 


303 


304 
305 


306 
307 
308 


E20 


Table E-l. (continued) 


TERR=IERR+2 
RETURN 
END 


SUBROUTINE LOWER(X,V,XL,XU,IERR) 
IMPLICIT REAL*8(A-H,0-Z) 

ae El -OF0T) GO 10 305 

TERR=0 

T=0 

TSET=X/V 

T=TSET 

TieOrs71SET 

Tee V.Gn.0. 14991 1=0.25*TSET 
TMAX=WOD( TSET) 


TST+TI 


[=I +] 

IF (1.£Q.20) GO TO 304 

IF (WOD(T)/TMAX.GE.0.001) GO TO 300 
XU=T 

T=0 

TL=TSET 

IF (V.GT.0.149) GO TO 303 

We lec 

T=I+] 

JEMti EO 20) GO 10) 306 

IF (WOD(TL)/TMAX.GE.0.001) GO TO 302 
XL=TL 

GO TO 307 

TL=0.75*TL 

T=I+] 

IF (1.EQ.20) GO TO 306 

IF (WOD(TL)/TMAX.GE.0.001) GO TO 303 
XL=TL 

GO TO 307 

NUS Ors Ell 

GO TO 301 

XU-1000 .0 

XL=X 

GO TO 308 

Ea) 20s Z071SET 

IERR=IERR+20 

IF (XU.GT.1000.0)XU=1000.0 

IF (XL.GT.1000.0)XL=1000.0 

RETURN 

END 
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Table E-1. (continued) 


FUNCTION FCT(T) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON ALPHA,BETA,AX,NX/FUN/VEL , XX 

POl=We0y iz*2. //0SORT (1. OFBETA*T**2) *DEXP( -(XX-VEL*T ) **2/ 
are Th O+BE TAs 2 ) ) ) 

RETURN 

END 


FUNCTION WOD(T) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON ALPHA,BETA,AX,NX/FUN/VEL , XX 

WOD=1 .0/T**1 .5/DSQRT(1 .0+BETA*T**2 ) *DEXP(-(XX-VEL*T)**2/ 
Ween )* 1 OtBeiAeie. 2) ) ) 

RETURN 

END 
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Table E-2. FORTRAN IV Computer Program for a 
Contained Reactor System 


AKEKKSE ADT FxxkK** 


CONTINUOUS RELEASE CASE 
WITH DECAY 


IMPLICIT REAL*8(A-H,0-Z) 

EXTERNAL FCT 

COMMON ALPHA,BETA,AX,NX/ISOT/IDENT(200) ,AMT(200) ,PCC(200) , 
IDIST(200) , TIME(200)/FUN/VEL ,XX/MISC/POWER, GAM(200, 7) ,PER(200) 
IRTOT(40) ,CONC(40) , VDIS(16) 


PRINT HEADING,READ AND PRINT POWER LEVEL 


PRINT 10 

10 FORMAT ('1*****SEADIF***** CONTINUOUS RELEASE CASE") 
READ 11,POWER 

11 FORMAT (F10.4) 
PRINT 12,PQWER 

12 FORMAT ('OPOWER=',F10.4,' MW') 


READ IN CURRENT INFO AND COMPUTE ALPHA AND BETA 


READ 13,AX,AY,AZ,OMEGAY ,OMEGAZ 
13 FORMAT (2F10.3,3D10.3) 
ALPHA=1 0/8.0/3.1415927**1 .5/DSORT (AX*AY*AZ) 
BETA=(OMEGAY **2*AY /AX+OMEGAZ**2*AZ/AX) 12.0 
PRINT 15,AX,AY,AZ,OMEGAY ,OMEGAZ,ALPHA, BETA 
15 FORMAT(' AX=',F10.3,' METERS SQUARED PER SEC'/ 
1' AY=',F10.5,' METERS SQUARED PER SEC'/ 
2' AZ=',F10.5,' METERS SQUARED PER SEC'/ 
3' OMEGA Y=',D10.3,' PER SEC'/ 
4' OMEGA Z=',D10.3,' PER SEC'/ 
5' ALPHA =',D10.3,' SEC*3/2 PER METER*3'/ 
6' BETA =',D10.3,' PER SEC SQUARED') 
READ 16,NV,NX 
16 FORMAT (213) 
READ 72 (VOU) 
READ 18,(x0(I) 
17 FORMAT (8F10.6 
18 FORMAT (8F10.2 
READ 19,RATE 
19 FORMAT (D10.3) 
PRINT 20,RATE 
20 FORMAT (' LEAKAGE RATE TO ENVIRONMENT=',D10.3,' PER SEC') 


»1=1,NV) 
,1=1,NX) 
) 
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Table E-2 (continued) 


READ IN EQUIVALENT CONTAINMENT DIAMETER (METERS) AND COMPUTE THE 
VIRTUAL DISTANCES FOR THE GIVEN VELOCITIES 


READ 21,DIAM 

FORMAT (F10.2) 

PRINT 22,DIAM 

FORMAT (' CONTAINMENT DIAMETER =',F10.2,' METERS') 
VIRT=DIAM**2/72 .0/AX**0. 75/AY**0 25 

DO 23 I=1,NV 

VDIS(I)=VO(I)*VIRT 


READ IN ISOTOPE INFO AND COMPUTE SOURCE TERM 
READ 24,NI 


FORMAT (13) 
READ 25,((IDENT(I),AMT(I) ,PCC(I) ,HALF(I),(GAM(I,J),J=1,7, 


IPER(I)),1=1,NI) 


FORMAT (15.5%, 3090.3 65.1 ,F5.3) 
CALL SOURCE(NI,XEM) 

PRINT 26,XEM 

FORMAT (' XEM='D10.3,' RM*3/SEC') 


COMPUTE DOSE RATE (STEADY STATE) 


PRINT 28 
FORMAT ('IRADIATION DOSAGE FOR STEADY STATE CONDITIONS'// 


1'  DIST(M) VEL(M/SEC) DOSE RATE(R/HR)  4HR DOSE(R)' 
2' INFINITE DOSE LIMIT VIRT DIST(M)'/) 


NINT=500 

DAX=DSQRT (AX) 

DO 47 I=1,NV 

VEL=VO(I) 

IF (VEL.EQ.0.0)NINT=1200 

DO 32 J=1,NX 

XX=X0(J)+VDIS(1) 

PF (VEEGE 051) GO 10) 31 

CALL LIMIT(XX,VEL,XL,XU,IERR) 
S=DGNTGR(XL,XU,NINT,FCT) 
TERM(I,J)=2.0*S 
RHR(J)=XEM*RATE*ALPHA*146.52*TERM( I,J) 
RTOT(J)=XEM*RATE*ALPHA*3. 066D+3*TERM( I J) 
RPOP (J )=XEM*RATE*ALPHA*6 .00D+3*TERM(I ,J) 
PRINT 30,X0(J),VEL,RHR(J) ,RTOT(J) ,RPOP(J) ,IERR, VDIS(1I) 
FORMAT (! SMO ac 3X, EOL 6 3 Shean Ss) 126456 10.5) 
GO TO 32 

S=3.5449*DAX/ XX 

IERR=0 

GO T0 29 

CONTINUE 
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Table E-2. (continued) 


DETERMINE EXCLUSION RADIUS 


IF (RTOT(1)-25.0) 34,34,36 
34 PRINT 35,V0(1),X0(1) 
35 FORVAiee GOR VELOCITY OF ',FI@26,' M/SEC, THE EXCLUSION' 
VRAEIUSSIS "LESS THAN®,FIO.3,° ‘METERS’ ) 


GO TO 40 

36 DO 38 K=2,NX 
IF (RTOT(K)-25.0) 37,38,38 

37 XCRE=(RTOT(K-1)-RTOT(K) )/(X0(K) -X0(K-1) ) 
PLUS=(RTOT(K-1)-25.0)/XCRE 
XCLU=X0(K-1)+PLUS 
PRINT 35,VO(I),XCLU 


GO TO 40 

38 CONTINUE 
PRINT 39,VO(I) ,XO(NX) 

39 RORVAT OFOR VELOCITY OF ’,FIOl6,° “M/SEC, TRE EXCLUSION’ 
1' RADIUS IS GREATER THAN',F10.3,' METERS') 


FIND LOW POPULATION RADIUS 


40 IF (RPOP(1).LT.25.0) GO TO 45 
DO 4] N=2,NX 
IF (RPOP(N).LT.25.0) GO TO 43 

41 CONTINUE 
PRINT 42,X0(NX) 

42 FORMAT(' ',35X,'THE LOW POPULATION RADIUS IS GREATER THAN’, 
1EIO. cee METERS.) 
GO TO 47 

43 XCRE=(RPOP(N- 
PLUS=(RPOP(N 
RAD=X0(N-1)+ 
PRINT 44,RAD 

44 FORMAT(' ',35X,'THE LOW POPULATION RADIUS IS APPROXIMATELY ', 
lEV@re same GERS:/) 
GO T0 47 

45 PRINT 46,X0(1) 

46 FORMAT (' ',35X,'THE LOW POPULATION RADIUS IS LESS THAN',F10.3,' 
IMETERS '/) 

47 CONTINUE 


1)-RPOP(N))/(XO(N)-XO(N-1) ) 
1)-25.0)/XCRE 
LU 


PLUS 


COMPUTE CONCENTRATIONS 


PRINT 48 

48 FORMAT ('ICRITICAL ISOTOPE INFORMATION'/' ',10X,'NOTE 1 -' 
1' INDICATES DISTANCE LESS THAN MINIMUM DISTANCE.'/11X, 
2'NOTE 2 - INDICATES DISTANCE GREATER THAN MAXIMUM DISTANCE. ' 
3//' ISOTOPE',3X,'CURIES/MW' 6X, 'MPCC(UC/KG) HALF LIFE(HRS)' 
4' — DISTANCE (M) NOTE'/) 
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Table E-2. (continued) 


DO 54 K=1,NI 

IF (PCC(K).EQ.0.0) GO TO 49 

DOTM=AMT(K) *POWER*PER(K) *RATE*ALPHA* 1000.0 

NOTE=0 

CONC(1)=DOTM*TERM(NV,1) 

EeveonerliiGE. PCC(K)) GO TO 50 

DIST(K)=X0(1) 

NOTE=1 

GO TO 52 

DIST(K)=0.0 

GO TO 54 

DO 51 N=1,NX 

CONC(N)= DOTM*TERM(NV ,N) *DEXP(- 0.693/HALF(K)*X0(N)/VO(NV) ) 

_ (CONC(N).GT.PCC(K)) GO TO 51 
F (N.EQ.1)DIST(K)=X0(1) 

iF (N.EQ.1) i [Ue Sy 

DIST(K)=XO(N-1)+(CONC(N~1)-PCC(K) )*(XO(N)-XO(N-1) )/(CONC(N-1)- 


1CONC(N) ) 


GOmi0 32 

CONTINUE 

DIST(K)=X0(NX) 

NOTE=2 

PRINT 53, IDENT(K) ,AMT(K) ,PCC(K) ,HALF(K) ,DIST(K) ,NOTE 
Bervede 515.3% .4(D11.3,5X), 12) 

CONTINUE 

STOP 

END 


SUBROUTINE SOURCE(NI,XEM) 
IMPLICIT REAL*8(A-H,0-Z) 
COMMON/ ISOT/IDENT(200) ,AMT( 200) ,PCC(200) ,DIST(200) , TIME (200) 


1/MISC/POWER ,GAM( 200, 7) ,PER(200) 


DIMENSION A(200,7),SUM(7) ,EN(7), TEMP(7) 
DO 100 I=1,NI 

DO 100 J=1,7 
A(1,J)=AMT(1)*GAM(1,J)*PER(1)*0.01*POWER 
CONTINUE 

DO 101 J=1,7 

SUM(J)=0.0 

DO 101 I=1,NI 

SUM(J)=SUM(J)+A(1,d) 


CONTINUE 

EN(1)=4.30D-15 
EN(2)=1.46D-14 
EN(3)=2.04D-14 
EN(4)=3.17D-14 
EN(5)=2.05D-14 
EN(6)=4.17D-14 
EN(7)=5.26D-14 
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Table E-2. (continued) 


XEM=0 .0 

DOM OZ let 5/ 

TEMP (J)=EN(J)*SUM(J) 
XEM=XEM+TEMP (J) 
CONTINUE 
XEM=XEM*3.7D+10 
RETURN 

END 


SUBROUTINE LIMIT(X,V,XL,XU,IERR) 

LIMIT ERROR MESSAGE HAS THE FOLLOWING MEANINGS: 

-O- UPPER AND LOWER LIMIT LESS THAN 4 HRS. NO ERROR 

-l- UNABLE TO ACHIEVE DESIRED ACCURACY ON UPPER LIMIT 

-2- UNABLE TO ACHIEVE DESIRED ACCURACY ON LOWER LIMIT 

-3- UNABLE TO ACHIEVE DEISRED ACCURACY ON EITHER UPPER 
OR LOWER LIMIT 


IMPLICIT REAL*8(A-H,0-Z) 
IERR=0 

T=Q 

TE Giaeo.0.0) GO 10 210 
Sela 

TMAX=FCT(TSET) 

IF (V.LT.0.060) GO TO 201 
IF (¥.LT.0.100) GO TO 202 
Pemee le. 0500 mcd eTOs203 
l= 052 1SEM 

GO TO 204 

XU=100.0*TSET 

IF (V.GE.0.060) XU=20.0*TSET 
IF (V.GE.0.100) XU=10.0*TSET 
IF (V.GT.0.500)XU=5.O*TSET 
IERR=1 

GO TO 206 


Se Uetesie «Eat Tf 


5. OaiSET 

GO TO 204 

WETS ll 

GO TO 204 

IAG Sg Sia] 

T=TSET 

T=T+TI 

T=[+] 

IF (1.EQ.20) GO TO 200 
IF (FCT(T)/TMAX.GE.0.001) GO TO 205 
XU=T 


a 


IF (V.GT.0.500) GO TO 208 
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Table E-2. (continued) 


SET LOWER LIMIT 


I=0 

TL=TSET 

TL=TL/2.0 

[=I+] 

IF (I.£Q.20) GO TO 211 
IF (FCT(TL)/TMAX.GE.0.001) TO TO 207 
XL=TL 

GO TO 210 

1-0 © 

TL=TSET 

TL=0.8*TL 

[=I+] 

IF (1.£0.20) GO TO 211 
IF (FCT(TL)/TMAX.GE.0.001) GO TO 209 
XL=TL 

ae Te 2 

XL=X 

XU=1. OD+6*X 

Garner 212 
XL=1.0/2.0**20*TSET 
IERR=IERR+2 

RETURN 

END 


FUNCTION FCT(T) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON ALPHA,BETA,AX ,NX/FUN/VEL, XX 

FGI =120)7 1 * Ino, USORTCT OBE AW 1" "2 )“DEXEi=( XX =uE* [yen2/ 


1(4.0*AX*T* (1 .0+BETA* I**2) ) ) 


RETURN 
END 
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Figure E-l. Data Card Layout for SEADIF 


//DATA. INPUT DD * 








One to five cards 





Omit card entirely for uncontained release 





| One card for each isotope 
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